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Abstract 



For perfect fluids with equation of state p = p{n,s), Brown Q gave an 
action principle depending only on their Lagrange coordinates a^{x) without 
Clebsch potentials. After a reformulation on arbitrary spacelike hypersurfaces 
in Minkowski spacetime, the Wigner-covariant rest-frame instant form of these 
perfect fluids is given. Their Hamiltonian invariant mass can be given in closed 
form for the dust and the photon gas. The action for the coupling to tetrad 
gravity is given. Dixon's multipoles for the perfect fluids are studied on the 
rest-frame Wigner hyperplane. It is also shown that the same formalism can 
be applied to non-dissipative relativistic elastic materials described in terms 
of Lagrangian coordinates. 
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Typeset using REVT^ 



I. INTRODUCTION 



Stability of stellar models for rotating stars, gravity-fluid models, neutron stars, accretion 
discs around compact objects, collapse of stars, merging of compact objects are only some 
of the many topics in astrophysics and cosmology in which relativistic hydrodynamics is the 
basic underlying theory. This theory is also needed in heavy-ions collisions. 

As shown in Ref. there are many ways to describe relativistic perfect fluids by means of 
action functionals both in special and general relativity. Usually, besides the thermodynam- 
ical variables n (particle number density), p (energy density), p (pressure), T (temperature), 
s (entropy per particle), which are spacetime scalar fields whose values represent measure- 
ments made in the rest frame of the fluid (Eulerian observers), one characterizes the fluid 
motion by its unit timelike 4-velocity vector field (see Appendix A for a review of the 
relations among the local thermodynamical variables and Appendix B for a review of covari- 
ant relativistic thermodynamics following Ref. |0]). However, these variables are constrained 
due to [we use a general relativistic notation: "; /x" denotes a covariant derivative] 

i) particle number conservation, (nf/^).^ = 0; 

ii) absence of entropy exchange between neighbouring flow lines, {nsU'^)-^ = 0; 

iii) the requirement that the fluid flow lines should be fixed on the boundary. 
Therefore one needs Lagrange multipliers to incorporate i) anf ii) into the action and this 
leads to use Clebsch (or velocity-potential) representations of the 4-velocity and action func- 
tionals depending on many redundant variables, generating first and second class constraints 
at the Hamiltonian level (see Appendix A). 

Following Ref. § the previous constraint iii) may be enforced by replacing the unit 4- 
velocity If^ with a set of spacetime scalar fields a^{z), i = 1,2, 3, interpreted as "Lagrangian 
(or comoving) coordinates for the fluid" labelling the fluid flow lines (physically determined 
by the average particle motions) passing through the points inside the boundary (on the 
boundary they are fixed: either the &{z°, z)^s have a compact boundary Va{z°) or they have 
assigned boundary conditions at spatial infinity). This requires the choice of an arbitrary 
spacelike hypersurface on which the a*'s are the 3-coordinates. A similar point of view is 
contained in the concept of "material space" of Refs. 0,^, describing the collection of all 
the idealized points of the material; besides to non-dissipative isentropic fluids the scheme 
can be applied to isotropic elastic media and anisotropic (crystalline) materials, namely to 
an arbitrary non-dissipative relativistic continuum ||^. See Ref. for the study of the 
transformation from Eulerian to Lagrangian coordinates (in the non-relativistic framework 
of the Euler-Newton equations). 

Notice that the use of Lagrangian (comoving) coordinates in place of Eulerian quantities 
allows the use of standard Poisson brackets in the Hamiltonian description, avoiding the 
formulation with Lie- Poisson brackets of Ref. |P, which could be recovered by a so-called 
Lagrangian to Eulerian map. 

Let be a curved globally hyperbolic spacetime [with signature e(H ), e = ±] 

whose points have locally coordinates z'^. Let '^g^u{z) be its 4- metric with determinant 
'^g = \det^g^y\. Given a perfect fluid with Lagrangian coordinates a{z) = {a*(z)}, unit 
4-velocity vector field U^{z) and particle number density n{z), let us introduce the number 
flux vector 
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n{z)U^'{z) = (1.1) 



and the densitized fluid number flux vector or material current [e*^^^^ = 1/a/^; 
da{\/^^^^^") = 0; rji2z{di') describes the orientation of the volume in thematerial space] 



"-(^J = = ^123 (tt U))- 



d,r{a\z)) = [n{z)U^{z)],^ = 0, 
^ J^(a*(;2))9^a^(^) = [V^nf/^](^)9^«^(;2) = 0. (1.2) 

This shows that the fluid flow lines, whose tangent vector field is the fluid 4- velocity timelike 
vector field U^, are identified by 5' = const, and that the particle number conservation is 
automatic. Moreover, if the entropy for particle is a function only of the fluid Lagrangian 
coordinates, s = s(a'), the assumed form of also implies automatically the absence of 
entropy exchange between neighbouring flow lines, {nsU^)^^ = 0. Since t/'^9^s(a'') = 0, the 
perfect fluid is locally adiabatic; instead for an isentropic fluid we have 9^5 = 0, namely 
s = const.. 

Even if in general the timelike vector field U^{z) is not surface forming (namely has 
a non- vanishing vorticity, see for instance Ref. |^), in each point z we can consider the 
spacelike hypersurface orthogonal to the fluid flow line in that point (namely we split the 
tangent space TjM*^ at z in the U^{z) direction and in the orthogonal complement) and 
consider ■^\U^e^i,p^dz^ A dz^ A dz°']{z) as the infinitesimal 3- volume on it at z. Then the 
3-form 

ri[z] = [r]i23{a)da^ A da^ A da^]{z) = ^n{z)[U^epupadz'' A dz'' A dz'']{z), (1.3) 

may be interpreted as the number of particles in this 3- volume. If V is a volume around z 
on the spacelike hypersurface, then Jy rj is the number of particle in V and /y sr] is the total 
entropy contained in the flow lines included in the volume V. Note that locally 77123 can be 
set to unity by an appropriate choice of coordinates. 
In Ref. [Ili it is shown that the action functional 



^[V, «] = - / d'z^^)p{^-^^^^M^, sia\z))), (1.4) 

has a variation with respect to the 4-metric, which gives rise to the correct stress tensor 
T^*^ = (p + p)U^U^ — ep^g^" with p = n^\s — p for a perfect fluid [see Appendix A]. 

The Euler-Lagrange equations associated to the variation of the Lagrangian coordinates 
are [Q [V^ = /it/^ is the Taub current, see Appendix A] 



/4 



^ ^^-d^a' = ^e"^^V„;^[/^e,^^5 - T d,s = 2%;,]f/'^ - T d,s = 0. (1.5) 



^/^Sa' " 2 
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As shown in Appendix A, these equations together with the entropy exchange constraint 
imply the Euler equations imphed from the conservation of the stress-energy-momentum 
tensor. 

Therefore, with this description the conservation laws are automatically satisfied and the 
Euler- Lagrange equations are equivalent to the Euler equations. In Minkowski spacetime 
the conserved particle number is A/" = Jva(^z°) d^zn{z) U°{z) = Jv^(^z°) ^^^^ while the 

conserved entropy per particle is Jy^f^^o-^ z s{z) n{z) U° (z) = Jy^^^^o^d^z s{z) J'^^a'-). More- 
over, the conservation laws T^^" ,^ = will generate the conserved 4- momentum and angular 
momentum of the fluid. 

However, in Ref. there are only some comments on the Hamiltonian description im- 
plied by this particular action. 

This description of perfect fluids fits naturally with parametrized Minkowski theories 



for arbitrary isolated relativistic systems [see Ref. |TI[] for a review] on arbitrary spacelike 
hypersurfaces, leaves of the foliation of Minkowski spacetime associated with one of its 
3+1 splittings. 

Therefore, the aim of this paper is to find the Wigner covariant rest-frame instant form 
of the dynamics of a perfect fluid, on the special Wigner hyperplanes orthogonal to the total 
4-momentum of the fluid. In this way we will get the description of the global rest frame 
of the fluid as a whole; instead, the 4-velocity vector field If^ defines the local rest frame 
in each point of the fluid by means of the projector '^g'^'^ — eU^lf^. This approach will also 
produce automatically the coupling of the fluid to ADM metric and tetrad gravity with the 
extra property of allowing a well defined deparametrization of the theory leading to the 
rest-frame instant form in Minkowski spacetime with Cartesian coordinates when we put 



equal to zero the Newton constant G |Tl|. In this paper we will consider the perfect fluid 
only in Minkowski spacetime, except for some comments on its coupling to gravity. 

The starting point is the foliation of Minkowski spacetime M^, which is defined by an 
embedding i? x E — >• M"^, (r, a) z'^{t,(t) G S,- and with S an abstract 3-surface dif- 
feomorphic to R^, with S,- its copy embedded in labelled by the value r (the scalar 
mathematical "time" parameter r labels the leaves of the foliation, a are curvilinear coor- 
dinates on and = (a'^ = T,a^) are ET--adapted holonomic coordinates for M"^). See 
Appendix C for the notations on spacelike hypersurfaces. 

In this way one gets a parametrized field theory with a covariant 3+1 splitting of 
Minkowski spacetime and already in a form suited to the transition to general relativity 



in its ADM canonical formulation (see also Ref. ||T2[ , where a theoretical study of this prob- 
lem is done in curved spacetimes). The price is that one has to add as new independent 
configuration variables the embedding coordinates z^{T,a) of the points of the spacelike 
hypersurface S,- [the only ones carrying Lorentz indices] and then to define the fields on E7- 
so that they know the hypersurface of r-simultaneity [for a Klein-Gordon field 0(x), this 
new field is 0(r, a) = 0(-2(t, a)): it contains the non-local information about the embedding]. 
Then one rewrites the Lagrangian of the given isolated system in the form required by the 
coupling to an external gravitational field, makes the previous 3+1 splitting of Minkowski 
spacetime and interpretes all the fields of the system as the new fields on (they are 
Lorentz scalars, having only surface indices). Instead of considering the 4- metric as de- 
scribing a gravitational field (and therefore as an independent field as it is done in metric 
gravity, where one adds the Hilbert action to the action for the matter fields), here one 
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replaces the 4-metric with the the induced metric QabI^] — z'^Vij.i'^'^ [a functional 

of z^; z'^ = dz'^/da^ are fiat inverse tetrad fields on Minkowski spacetime with the z^'s 
tangent to S,-] and considers the embedding coordinates z^{T,a) as independent fields [this 
is not possible in metric gravity, because in curved spacetimes z'^ ^ dz^/ da^ are not tetrad 
fields so that holonomic coordinates z^{t, a) do not exist]. Prom this Lagrangian, besides a 
Lorentz-scalar form of the constraints of the given system, we get four extra primary first 
class constraints 

n.ir, a) = p,{t, a) - 1,{t, a)Tj;,{r, a) - z,,{t, a)Tg{r, a) ^ 0, (1.6) 

[ here Tgyg{T,a) = A4{T,a), T^yg{T,a) = A4^{T,a), are the components of the energy- 
momentum tensor in the holonomic coordinate system, corresponding to the energy- and 
momentum-density of the isolated system; one has {H(fM){T,a),Ti.(u)iT,ff)} — 0] implying 
the independence of the description from the choice of the 3+1 splitting, i.e. from the choice 
of the foliation with spacclike hypersufaces. As shown in Appendix C the evolution vector 
is given by 2;^ = N[z]{fiat)^^ + ^[z]{}iat)^f 1 where l^{r, a) is the normal to S^- in z^{t, a) and 
N[z]{fiat){Ti (^), ^[z\{fiat){^^ ^) lapsc and shift functions defined through the metric 

like in general relativity: however, now they are not independent variables but functionals 
of z^'{T, a). 

The Dirac Hamiltonian contains the piece / dPa\^{T, a)l-L^(T^ a) with \^{t, a) Dirac mul- 
tipliers. It is possible to rewrite the integrand in the form [7'''' = —e^g^^ is the inverse of 
the spatial metric gfs = '^Qf-s = —^^drs, with ^^f^^^ of positive signature (-1- -|- -|-)] 

= Nuiat){T, o){l,n^){T, a) - N^fiat)f{r, a){Y'zs.nn{r, a), (1.7) 

with the (non-holonomic form of the) constraints (/^7i'*)(r, a) f» 0, {Y^Zs^Ti.'^){T,a) ^ 0, 
satisfying the universal Dirac algebra of the ADM constraints. In this way we have defined 
new fiat lapse and shift functions 

N^fiat){r,a)^X,{T,a)nT,a), 

N{fiat)f{T, B) = A^(t, B)z^{t, a). (1.8) 

which have the same content of the arbitrary Dirac multipliers Aj^(t, a), namely they multiply 
primary first class constraints satisfying the Dirac algebra. In Minkowski spacetime they are 
quite distinct from the previous lapse and shift functions N^z\{fiat), N[z]{fiat)f, defined starting 
from the metric. Since the Hamilton equations imply z!^{t, a) = X^{t, a), it is only through 
the equations of motion that the two types of functions are identified. Instead in general 
relativity the lapse and shift functions defined starting from the 4-metric are the coefficients 
(in the canonical part He of the Hamiltonian) of secondary first class constraints satisfying 
the Dirac algebra independently from the equations of motion. 

For the relativistic perfect fluid with equation of state p = p{n, s) in Minkowski space- 
time, we have only to replace the external 4-metric '^Qf^,, with g^^{T,a) = '^QAsi'^^^) ^i^d 
the scalar fields for the Lagrangian coordinates with Q;*(r, a) = a\z{T,a))] now either the 
Q;'(r, (t)'s have a compact boundary 14 (r) C or have boundary conditions at spatial in- 
finity. For each value of r, one could invert a* = Q!*(t, a) io a = (t{t, a*) and use the «*'s as 
a special coordinate system on inside the support Va{r) C S,-: z'^{t, a{T, a*)) — z^{t, a*). 
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By going to S^-adapted coordinates such that 77123 («) = 1 we get [7 = \det gfs\', y/g 



V^=^\detg^^\ = N^] 

j\a\T,a)) = [N^nU\r,a), 

r{a\T,a)) = [-e™*9,a^afia'9ca3](r,a) = -det {dfa'){r, a), 
3 

i=l;i,j,k cyclic 



with M = Jy^^^^d^aJ^{a'{T,a)) giving 

the conserved particle number and Jy^^^^-^d^als J'^){T,a) giving the conserved entropy per 
particle. 

The action becomes 

S = J dTd^aL{z^'{T,a),a\T, a)) = 



- 1 dTd^a{N^){T,a) 



1 



pi I \ iJ^)^-^9uv 77 77 iT,a;a'{T,a)), s{a\T,a))) 



N N 



;i.io) 



with N = N[,]^j, at), = N(^^^f^^,y 

This is the form of the action whose Hamiltonian formulation will be studied in this 
paper. 

We shall begin in Section II with the simple case of dust, whose equation of state is 
p = fin. 

In Section III we will define the "external" and "internal" centers of mass of the dust. 

In Section IV we will study Dixon's multipoles of a perfect fluid on the Wigner hyperplane 
in Minkowski spacetime using the dust as an example. 

Then in Section V we will consider some equations of state for isentropic fluids and we 
will make some comments on non-isentropic fluids. 

In Section VI we will define the coupling to ADM metric and tetrad gravity. 

In Section VII we will describe with the same technology isentropic elastic media. 

In the Conclusions, after some general remarks, we will delineate the treatment of perfect 
fluids in tetrad gravity (this will be the subject of a future paper). 

In Appendix A there is a review of some of the results of Ref. for relativistic perfect 
fluids. 
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In Appendix B there is a review of covariant relativistic thermodynamics of equihbrium 
and non-equihbrium. 

In Appendix C there is some notation on spacelike hypersurfaces. 

In Appendix D there is the definition of other types of Dixon's multipoles. 
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II. DUST. 



Let us consider first the simplest case of an isentropic perfect fluid, a dust with p = 0, 
s = const., and equation of state p = fin. In this case the chemical potential p is the rest 
mass-energy of a fluid particle: p = m (see Appendix A). 

Eq.( |l.l(J| ) implies that the Lagrangian density is [we shall use the notation gj^^ = ^g^B 

with signature e(H ), e = ±1; by using the notation with lapse and shift functions 

given in Appendix C we get: grr = €{N^ — ^grsN^N"^), grf = —e^grsN"^, gfs = —e^gfs with 
^gfs of positive signature (+++), g^^ = 5'"' = -e^, 9^' = -e{^9^' - ^jff^)] the inverse 
of the spatial 4-metric ^gfs is denoted 7*"* = ^7^'^ = —e^g'^'^, where ^g^'^ is the inverse of the 
3-metric ^gfs and we use ^7 = y/^gfs\ 



pNMry - ^grsY'Y^ = -pNX, 



X = ^{J^? - ^grsY'Y'^ = = V7^' (2.1) 



with J^, given in Eqs.(1.9). 



The momentum conjugate to a* is 

= /i— ^gif e™eyfc dua^d^a'' = P^ Tfi, 

where adJif = {detJ).!^^^ is the adjoint matrix of the Jacobian J = (Jif = dfo") of the 
transformation from the Lagrangian coordinates q;*(t, a) to the Eulerian ones a on S,-. 
The momentum conjugate to is 

Pf^i^,^) = ~^(^'^) = [P^f^^^ + P^r,,J^ j^\{r,a). (2.3) 
The following Poisson brackets are assumed 

K(r,a),p,(r,a'} = -<5=^(a-a'), 

{a\T,a),Uj{T,a')}=6';6%a-a'). (2.4) 
We can express Y"^ /X in terms of Hj with the help of the inverse (y-i)" of the matrix 

Tu 

^ 1 (2.5) 

A /i 
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where 



det {d^a^) 

Prom the definition of X we find 



iT-T - (2.6) 



X = ^-^^ (2.7) 

//.2 + 3^..(T-i)-(T-i)^m,n,- 



Consequently, we can get the expression of the velocities of the Lagrangian coordinates in 
terms of the momenta 

namely 

dra' = dfa' [n' - N{T-y'Iii ^/.2 + 3^..(T-i)-(T-i)'^m,n,] . (2.9) 
Now Pfj_ can be expressed as a function of the z's, a's and H's: 

Pm = J' V/^2 + 3to(T-i)-n,(T-i)*m,- + (T-^)"n,. (2.10) 

Since the Lagrangian is homogenous in the velocities, the Hamiltonian is only 

Hd = j d^aX^ir, a)n^{r, a), (2.11) 

where the Tif^ are the primary constraints 

^1, = Pi,-li,M + Zf^M^ ^ 0, 



satisfying 



M = T^^ = r ^fi^ + 3(^^i,(T-i)«n,(T-i)«m,-, 

= T^' = r {T-^Y'Iii. (2.12) 



{n^{T,a),n,{T,a)} = Q. (2.13) 



One finds that {7Y^(r, ci), H^} = 0. Therefore, there are only the four first class constraints 
?i^(r, a) ~ 0. They describe the arbitrariness of the foliation: physical results do not depend 
on its choice. 

The conserved Poincare generators are (the suffix "s" denotes the hypersurface E,-) 
= / d'ap''{r,a), 

Jr = / rfV[z^(r, a)p^(r, a) - z'^ir, a)p^{T, a)], (2.14) 

and one has 



K(r,a),J9r} = -r/^^ (2.15) 



[r, a) + 



+ / (fa 



Zf,J^{T'^rii,\{T,a)^Q. (2.16) 

Let us now restrict ourselves to spacelike hyperplanes S,- by imposing the gauge-fixings 
C'^(r, a) = z^{t, a) - x^(r) - 6^(r)a^ ^ 0, 

a),nAr, a)} = -r^^J\a - a'), (2.17) 

where x^(r) is an arbitrary point of S,-, chosen as origin of the coordinates a^, and bfi^r), 
f = 1,2, 3, are three orthonormal vectors such that the constant (future pointing) normal 
to the hyperplane is 

lf^(r, a)^F = K = e\p^h'({T)hl{T)hl{T). (2.18) 

Therefore, we get 

z>^{r,a)^h^{r), 
<(r,a)^x^(r)+6^(r)a^ 

grs{r, (?) ^ -t5rs, i'\t, a) ^ -e6'', 7(r, a) ^ 1. (2.19) 
By introducing the Dirac brackets for the resulting second class constraints 

{A, BY = {A, Bj-J d'a[{A, C^(r, a)}{n,iT, a),B}- {A, n,{r, a)}{e(r, a), B}], 

(2.20) 

we find that, by using Eq.(|2J[5|) and (|2l6| ) [with x^(r) = z'^(r, a) - b'^{T)a^ - (''{t, a) and 
with the assumption {6^(r),]9g} = 0], we get 

K(r),pr(r)}* = -v"". (2.21) 

The ten degrees of freedom describing the hyperplane are a;^(r) with conjugate momen- 
tum and six variables 0a (t), A = 1, .., 6, which parametrize the orthonormal tetrad &^(t), 
with their conjugate momenta Tx{t). 

The preservation of the gauge- fixings C^ir, a) ^ in time implies 

^Cir, ff) = mr, a), Hd} = -A^(r, a) - x^;(r) - b^{ry ^ 0, (2.22) 

so that one has [by using bi^: = and bf{T)b^{r) = — 6^(r)6^(r)] 

A^(r,a)^A'^(r)+A^(r)6,qr)(7^ 
A^(r) = -i;^(r), 

A'^^(r) = -r^(r) = Imrmr) - b'^irMr)]. (2.23) 
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Thus, the Dirac Hamiltonian becomes 



(2.24) 



and this shows that the gauge fixings C^{''~j (?) ~ do not transform completely the con- 
straints 'H^(r, (t) ^ in their second class partners; still the following ten first class con- 
straints are left 



- mT)K-b^,{T)b^,] J d'aa'{r^f,^ + 6^,{T-^rU,{T-^rm^){r,a) + 
+ [b'^imir) - b';{r)bUr)] J d'aa' {.r iT-y%)ir, a) ^ 0. (2.25) 
Here 5'^'^ is the spin part of the Lorentz generators 

•^S ~ -^sPs ~ -^sPs ~^ 5 

= 6^'(r) J d^aa'p-'ir, a) - b'^^r) J d^aa^'p^'ir, a). (2.26) 



As shown in Ref. ||T0[ instead of finding (j)x{T),Tx{T), one can use the redundant variables 



6^(r), S^'^{t), with the following Dirac brackets assuring the validity of the orthonormality 

condition r]''^ - b'^rj^h'^g = [C:^^^ = V^VfV^'^ + vf^VsV"" - VjVsV'"' - V^^Vfv'^'^ are the 
structure constants of the Lorentz group] 



(2.27) 



so that, while T-C^{t) ~ has zero Dirac bracket with itself and with T-C^^^t) ^ 0, these last 
six constraints have the Dirac brackets 



(2.28) 



We have now only the variables: x^, p^, 6^, S'^'^ , a*, Ilj with the following Dirac brackets: 

K(r),pr(r)r = -r/'^^ 

{5r(r),6^(r)}* = ryP-6^(r) - r/^^6;^(r), 
{5r(r),5f(r)r = C^7"^5f(r), 



K(r,a),n,(r,a)r = 5;.5^(a-a 



(2.29) 



After the restriction to spacelike hyperplanes we have 2^(r, a) ^ bf{T), so that zi^{T, a) ^ 
A^W(//at)(r,a)P(r,<T) + Ar[^j(^,^,)(r,cr) blr\T,a) ^ x^(r) + 6^(r)a^ = -A^(r) - >-(r)6,,(r)a^ 
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As said in the Introduction only now we get the coincidence of the two definitions of fiat 
lapse and shift functions (this point was missed in the older treatments of parametrized 
Minkowski theories): 

N[,]ifiat){r, a) ^ N^fiat){T, a) = -A^(r)P - P = N{t, a), 

(2.30) 

Let us now restrict ourselves to configurations with ep^ > and let us use the Wigner 
boost L'^i,(ps,Ps) to boost to rest the variables 6^, S^" of the following non-Darboux basis 



x^,P^,b%S^'^,a\Ili 



of the Dirac brackets {■,■}*■ The following new non-Darboux basis is obtained [x^ is no 
more a fourvector; we choose the sign rj — signp° positive] 



dp, 



SIX 



X 



M _ 



1 _ '^^''n'' 

Ps + \^Ps 



+ 



Ps + 



ltpl{p° 



-stp:i 



P's = 



a = a , 



= - le^{u{ps))VAB 



= -5r + 



1 



'^pUps + \/^pI 

We have 



dps, dp,^ ^^^^^ - 

■MsI^p: - stp':) + xf^siSTp: - sTp's)], 



(2.31) 



{^^prr = 0, 



Ps + y^Pl 



(2.32) 
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and we can define 

= ei{u{p,))e^{u{p.,))S^^ ^ [bf{r)b? - bf{r)b^] 

- [bHr)bf{r)-bfir)bfir)] J d'a a\riT-y%)ir, a). (2.33) 

Let us now add six more gauge-fixings by selecting the special family of spacelike liyper- 
planes S^vf ortliogonal to (tliis is possible for ep^ > 0), which can be called the 'Wigner 
foliation' of Minkowski spacetime. This can be done by requiring (only six conditions are 
independent) 

= blir) - eL^(«(Ps)) ^ 

bi{r)=ei{u{psW^{T)^r^i. (2.34) 

Now the inverse tetrad 6^ is equal to the polarization vectors e'X{u{ps)) [see Appendix 
C] and the indices 'f' are forced to coincide with the Wigner spin-1 indices 'r', while o = r 
is a Lorentz-scalar index. One has 

J d'a(.r K (T-YUi-a' iT-'Y%])ir, a), 

S7 ~ -S7 = - j d'cr{r a'- v'/i2 + <5„,(r-i)»n,(T-i)-^n,) (r, a). (2.35) 

The comparison of S^^ with yields 

guv ^ ^ur^vtgrt 

Xvr art f 

n I / 9 

The time constancy of ~ with respect to the Dirac Hamiltonian of Eq. (|2.24|) gives 

j^mr) - etf(«(p.))] = mr) - e'^iuips)), Hn}* = 

= lx-^iT){b^iT),Ssapir)r = A'^"(r)6,„(r) ^ 

A^''(r) ^ 0, (2.37) 

so that the independent gauge- fixings contained in Eqs.( |2.3^ ) and the constraints 'H'^^{t) ^ 
form six pairs of second class constraints. 

Besides Eqs.( p.l9| ), now we have [remember that x'^{t) = — A^(r)] 
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z'^ir) = i^(r) = 49i^)u^{ps) - Xs.{T)e^MPsMHPs)\ 

N{t) = ^) = [xs^{t)u^{Ps))], Vt = 1, 

_-2 __3 __r 

fl'ri -^s) 9rs — ^ fl'rs — ^^rsj 



[Xs ■ u{ps)Y 



(2.38) 



On the hyperplane all the degrees of freedom z'^{t, a) are reduced to the four degrees 
of freedom 5^ (r), which replace x^. The Dirac Hamiltonian is now = \^{t)'H^{t) with 



- j d^a{r[ui'{ps)\liJ? + 5„,(r-i)-n,(r-i)-m,- - 



0. 



(2.39) 



To find the new Dirac brackets, one needs to evaluate the matrix of the old Dirac brackets 
of the second class constraints (without extracting the independent ones) 



C 



^5^B[r^'^Ce%{u{ps))-v''''e^BHps))] 

{Tj, w^Y = {Tj, r|}* = 

\ = SAA[v'^e'A{u{Ps))-v''el{u{ps))] 



(2.40) 



Since the constraints are redundant, this matrix has the following left and right null 



eigenvectors: 



— 0,.3a 









[a^/j arbitrary], I ^, , j. Therefore, one has to find a left 

and right quasi-inverse C, CC — CC — D, such that C and D have the same left and right 
null eigenvectors. One finds 



C 



[Va.e^{u{ps)) - Vai^e^Hps))] 0^^ 



and the new Dirac brackets are 



n.al3D 

u{ps))eAT{u{ps)) 



(2.41) 



{^,5}- = {A,BY - \[{A,W'Y[ri,re^{<Ps)) - mre^HpM{Th.BY + 
+ {A,T^}*[r^,,eJ(«(p,)) -r^.^ef («(p,))]{71:'^^i?}*]. 



(2.42) 



While the check of {TY^^, BY* = is immediate, we must use the relation b^^T^e^P = -Tj 
[at this level we have = T^] to check {T^, BY* = 0. 
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Then, we find the following brackets for the remaining variables x^,p^, a*, Ilj 



and the following form of the generators of the "external" Poincare group 



(2.43) 



sJ^s 
oi 



Ps + \lm 

^iv ^js grs 



(2.44) 



Let us come back to the four first class constraints 1-L^{t) ^ 0, {Ti^ — O5 of 
Eg. ( p. 251) . They can be rewritten in the following form [from Eqs.( |1.9|) , ( p.6|) we have 



sys 



0, 



M. 



sys 



d^aM{T, a) 



J d'a{r^fi^ + (5„,(T-i)»n,(T-i)-^n,)(r, a 
- I d^a \det [drO^ 



\ 



"^pl^) Psys 



d^aM'\T, a) = 
d'afi(r{T-'Y%)ir, a) 



[det {dra 



d^a /i 



[r, a) ^ 0, 



(2.45) 



where Msys is the invariant mass of the fluid. The first one gives the mass spectrum of the 
isolated system, while the other three say that the total 3-momentum of the N particles on 
the hyperplane T^tW vanishes. 

There is no more a restriction on in this special gauge, because u^{ps) = 
gives the orientation of the Wigner hyperplanes containing the isolated system with respect 
to an arbitrary given external observer. Now the lapse and shift functions are 



-A(r) = a;^(r)M^(p,), 



(2.46) 



(2.47) 



so that the velocity of the origin of the coordinates on the Wigner hyperplane is 

x^(r) = e[-A(r)n^(p,) + A,(r)e^(«(p,)), [u'{p,) = 6, e'MPs)) = -e]. 
The Dirac Hamiltonian is now 

Hd = Hr)n{r) - A(r) • 7?p(r), (2.48) 
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and we have = {x^, H^}** = — A(r)'u''(ps)- Therefore, while the old had a velocity x^ 
not parallel to the normal = u^{ps) to the hyperplane as shown by Eqs. (|2.47|) , the new 
x^ has Xg ||/^ and no classical zitterbewegung. Moreover, we have that = / ■ x^ = / ■ x^ is 
the Lorentz-invariant rest frame time. 

The canonical variables x^, p^, may be replaced by the canonical pairs = -y/pf, Tg = 

Ps ■ isl^s [to be gauge fixed with - r ^ 0]; = Ps/e^ = u{ps), Zs = (^s{xs - ?fx°) = e^g;. 

One obtains in this way a new kind of instant form of the dynamics, the "Wigner- 
covariant 1-time rest-frame instant form" with a universal breaking of Lorentz covariance. 
It is the special relativistic generalization of the non-relativistic separation of the center of 
mass from the relative motion [H = ^ + Hrei] ■ The role of the "external" center of mass 
is taken by the Wigner hyperplane, identified by the point x^(r) and by its normal p^. The 
invariant mass M^ys of the system replaces the non-relativistic Hamiltonian Hrei for the 
relative degrees of freedom, after the addition of the gauge-fixing — r ^ [identifying 
the time parameter r, labelling the leaves of the foliation, with the Lorentz scalar time of 
the "external" center of mass in the rest frame, Tg = Ps ■ Xs/Msys and implying A(r) = — e]. 
After this gauge fixing the Dirac Hamiltonian would be pure gauge: Hr, = — A(r) • Hp{r). 
However, if we wish to reintroduce the evolution in the time r = in this frozen phase 
space we must use the Hamiltonian [in it the time evolution is generated by Mgys'- it is like 
in the frozen Hamilton- Jacobi theory, in which the evolution can be reintroduced by using 
the energy generator of the Poincare group as Hamiltonian] 

HD = Msys-X{T)-n,{T). (2.49) 

The Hamilton equations for a*(r, a) in the Wigner covariant rest-frame instant form are 
equivalent to the hydrodynamical Euler equations: 



dra\T,a) = {a\T,a),HD} = 



+ X'iT)dra\T), 



r, a) + 



drlliir, a) 



{n, 

_d_ 

da' 



r, cr 



{ijk cyclic) 



[det {dra'')] 



+ 



,,2 I ^uv duQ-"^dva^ tt tt 



[r, a) + 



_d_ 

da' 



det {dta^ 



-n„, + 



d a™ 



(r,a). 



(2.50) 



In this special gauge we have 6^ = L^AiPs,Ps) (^he standard Wigner boost for timelike 
Poincare orbits), S^" = S^^, [S^y, = e™'' / c/V a" (j^ (T-i)^^H,) (r, ct)], and the only re- 
maining canonical variables are the non-covariant Newton- Wigner-like canonical "external" 
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— * 

3-center-of-mass coordinate Zg (living on the Wigner hyperplanes) and kg. Now 3 degrees of 
freedom of the isolated system [an "internal" center-of-mass 3- variable g^ys defined inside the 
Wigner hyperplane and conjugate to Psys\ become gauge variables [the natural gauge fixing 
to the rest-frame condition P^yg ^ is X^yg ^ 0, implying \r[j) = 0, so that it coincides 
with the origin a;^(r) = z'^{t, a = 0) of the Wigner hyperplane]. The variable x'^ is playing 
the role of a kinematical "external" center of mass for the isolated system and may be inter- 
preted as a decoupled observer with his parametrized clock (point particle clock). All the 
fields living on the Wigner hyperplane are now either Lorentz scalar or with their 3-indices 
transformaing under Wigner rotations (induced by Lorentz transformations in Minkowski 
spacetime) as any Wigner spin 1 index. 
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III. EXTERNAL AND INTERNAL CANONICAL CENTER OF MASS, 
MOLLER'S CENTER OF ENERGY AND FOKKER-PRYCE CENTER OF 

INERTIA 



Let us now consider the problem of the definition of the relativistic center of mass of 
a perfect fluid conflguration, using the dust as an example. Let us remark that in the 
approach leading to the rest-frame instant form of dynamics on Wigner's hyperplanes there 
is a splitting of this concept in an "external" and an "internal" one. One can either look 
at the isolated system from an arbitrary Lorentz frame or put himself inside the Wigner 
hyperplane. 

From outside one flnds after the canonical reduction to Wigner hyperplane that there 
is an origin a;^(r) for these hyperplanes (a covariant non-canonical centroid) and a non- 
covariant canonical coordinate x^(r) describing an "external" decoupled point particle ob- 
server with a clock measuring the rest-frame time Tg. Associated with them there is the 
"external" realization ( p.44[ ) of the Poincare group. 

Instead, all the degrees of freedom of the isolated system (here the perfect fluid conflg- 
uration) are described by canonical variables on the Wigner hyperplane restricted by the 
rest-frame condition Pgyg ~ 0, implying that an "internal" collective variable (jsys is a gauge 
variable and that only relative variables are physical degrees of freedom (a form of weak 
Mach principle). 

Inside the Wigner hyperplane at r = there is another realization of the Poincare 
group, the "internal" Poincare group. Its generators are built by using the invariant mass 
Msys and the 3- momentum Pgys, determined by the constraints ( [^.251 ), as the generators of 
the translations and by using the spin tensor S^^ as the generator of the Lorentz subalgebra 



pr 
pr 

r 



M. 



sys 



P. 



sys 
jrr __ 



r,(T 



a 



S7^ 



S\ 



sys 



2 



r, o- 



0, 



(3.1) 



By using the methods of Ref. [0 (where there is a complete discussion of many deflni- 
tions of relativistic center-of-mass-like variables) we can build the three "internal" (that is 
inside the Wigner hyperplane) Wigner 3-vectors corresponding to the 3-vectors 'canonical 
center of mass' qsys, 'MoUer center of energy' fsys and 'Fokker-Pryce center of inertia' ygyg 
[the analogous concepts for the Klein-Gordon fleld are in Ref. |T^ (based on Refs. [0), 
while for the relativistic N-body problem see Ref. fl^ and for the system of N charged 
scalar particles plus the electromagnetic fleld Ref. p^ j. 

The non-canonical "internal" M0ller 3-center of energy and the associated spin 3-vector 

are 



sys 



K 
p7 



2P^ 



a a 



JV/i^ + 5™(T-i)«(T-i)''in,n, 



r, (T 



18 



^sys ^ '^sys ^ -^7 

pr 

{.^sysy } ^ ' {.^sys^ } p-j- ' 

V syst ' sysf ^pT^2 sysi 

{^'sys, ^lys} = ^''\Ks ' J^^i^^ys ' P) P% {^lys. = 0- (3-2) 

The canonical "internal" 3-center of mass qsys [{(fsys-, <fsys} = 0; {llys-iP^} = ^^^i 



Q.sys ^sys 



J X Vlgyg 



X J X P 

+ , , + 



^(pr)2 _ p2 ^(pr)2 _ p2(pr + ^(pr)2 _ ^2) 

(A' • P) P 

pr^(pr)2_p2(^pr + ^(pr)2 _ p2^ ' 

— * 

~ ^sj/s /or P f5i 0; {?sj/s, -P^} = ^ ~ 0, 



sys ^ Qsys ^ 



P^ J ^ K X P (J-P) P 



^{Pry - p2 ^(pr)2 _ p2 y'(pr)2 _ p2 (^pr + ^(pr)2 _ p2^ 

«4„ /or P«0, 5r,, = e-"|(^Va"(j-(T-irn,)(r,a), 

{'S'gsj/s, P} = {'S'qrsj^s, gsj^s} = 0, {'S'gsj^s, 'S'^sj^s) = ^^^^Sqsys- (3-3) 

The "internal" non-canonical Fokker-Pryce 3-center of inertia' y^j^s is 

Sgys X P Sgys X P 

Z/sj/s — Qsys H I ^ I fsys H y i 



(pr)2 _ p2(pr + J(pr)2 _ p2) P^ ^ {P^f - P^ 



Ssys X P PV,,, + V {P^y - P^Vsys 

Qsys ^sys + / — ~ I — i 

p^{p^ + \l{P^f - p2) p^ + v(P^)2 - p2 



\y sysi iysysS j sys ^ 



{Sgys ' P) P^ 



pr^(^prY_p2 - ^{P-y-P\P^ + ^{P-y-P^y 

P fs =^ Qsys ~ r^ys ~ Usys- (3-4) 

The Wigner 3- vector Qgyg is therefore the canonical 3-center of mass of the perfect fluid 
configuration [since Qsys ~ Tgyg, it also describe that point z'^{t, qsys) — 3:^(t) -|- qlyge!;!:{u(ps)) 
where the energy of the configuration is concentrated]. 
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There should exist a canonical transformation from the canonical basis a*(T, a), Ili{T, a), 
to a new basis Qsys, P = P^, al^iij^a), nj,eZi(r, a) containing relative variables a*g;(r, a), 
nre/i(T, c?) with respcct to the true center of mass of the perfect fluid configuration. To 



identify this final canonical basis one shall need the methods of Ref. ||16|| . 

The gauge fixing q^ys ~ [it implies A(r) = 0] forces all three internal center-of- 
mass variables to coincide with the origin of the Wigner hyperplane. We shall denote 
= x^{Qi) + TU^{ps) the origin in this gauge (it is a special centroid among the many 
possible ones; a;^(0) is arbitrary). 

As we shall see in the next Section, by adding the gauge fixings Xgy^ = qsys ~ one 
can show that the origin x^(r) becomes simultaneously the Dixon center of mass of an 



extended object and both the Pirani and Tulczyjew centroids (see Ref. for a review of 
these concepts in relation with the Papapetrou-Dixon-Souriau pole-dipole approximation of 
an extended body). The worldline x^^"'""'^^ is the unique center-of-mass worldline of special 



relativity in the sense of Refs. [19 



With similar methods from the rest-frame instant form "external" realization of the 
Poincare algebra of Eq. (|2.44|) with the generators p^, J^^ = x^pi — x{p\ + Kl = 

J°* = x°pI — xIp° = x"pI — xIp° + pt^^e [^°^ x° = this is the Newton- Wigner 

decomposition of J^'^] we can build three "external" collective 3-positions (all located on 
the Wigner hyperplane): i) the "external canonical 3-center of mass Qs connected with the 
"external" canonical non-covariant center of mass x^; ii) the "external" M0ller 3-center of 
energy Rg connected with the "external" non-canonical and non-covariant M0ller center of 
energy i?^; iii) the "external" Fokker-Pryce 3-center of inertia connected with the "external" 
covariant non-canonical Fokker-Price center of inertia y/* (when there are the gauge fixings 
ffsys ~ it coincides with the origin x^). It turns out that the Wigner hyperplane is the 



natural setting for the study of the Dixon multipoles of extended relativistic systems [|20 
(see next Section) and for defining the canonical relative variables with respect to the center 
of mass. 

The three "external" 3- variables, the canonical Qs, the M0ller Rg and the Fokker-Pryce 
Yg built by using the rest-frame "external" realization of the Poincare algebra are 

Rg = --Kg = (I - ^xl^ 



f) _H_^™0_^_pi "^sys 'X Ps _ P"Ps + ^sYs 

Xg n S -5 ' 



P°s ' ^s P'siPt + ^s) P^s + eg 

Ssys ^ Ps ^sys ^ Ps 



Yg = Qg + ^^^^ = Rs + 



{R:, Rl} = --i-e--^:, ng = Jg-RgX pg 

{PsY 

■ Ps) 

l-'s5-'sJ ,„n '-'sys ' 



^sPg 



^s{p°g + eg) 



Ps ' Qs Ps ' -^s Ps ' '^s ^s ' ^s 
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= ^ 4 = n = Rs, (3.5) 



with the same velocity and coinciding in the Lorentz rest frame where = €^(1; 0) 

In Ref. fl^ in a one-time framework without constraints and at a fixed time, it is shown 
that the 3- vector Ys [but not Qg and Rg] satisfies the condition {Kl,Yf} = YJ {Yf,p°} for 
being the space component of a 4-vector Y"/*. In the enlarged canonical treatment including 
time variables, it is not clear which are the time components to be added to Qg, Rs, Yg, to 
rebuild 4-dimesnional quantities x^, i?^, Y"/*, in an arbitrary Lorentz frame F, in which the 
origin of the Wigner hyperplane is the 4-vector = {x°; Xg)- We have 



1 r „„„ ^^..PsuP'^. 



(x^ir); x,{r)) = - . , ^ , PsuS^ + esiST - S] 

^S [Ps ' ^ 

kf. 



s ^2 



Pi 



= Vl + kliT, + ^-^) = Vl + kliT, + K ■ qs) + = 6,Vl + ^1 



es 

Xs = — + (Ts + ^^-^)ks = qs + (Ts + ks ■ qs)ks, Ps = esh. (3.6) 

for the non-covariant (frame-dependent) canonical center of mass and its conjugate momen- 
tum. 

Each Wigner hyperplane intersects the worldline of the arbitrary origin 4-vector x^(r) = 
z'^{t, 0) in (T = 0, the pseudo worldline of a;^(r) = z'^{t, a) in some a and the worldline of 
the Fokker-Pryce 4-vector (t) = z^{t, ay) in some ay [on this worldline one can put the 
"internal center of mass" with the gauge fixing ~ (g^ ^ ^ due to ~ 0)]; one 
also has i?^ = z'^{t,(Tr). Since we have Tg = u{ps) ■ Xg = u{ps) ■ Xg = t on the Wigner 
hyperplane labelled by r, we require that also Y"/*, i?^ have time components such that they 
too satisfy u{ps) ■ Yg = u{ps) ■ Rg = = t. Therefore, it is reasonable to assume that x^, 
Y^ and i?^ satisfy the following equations consistently with Eqs.(|3.2|), ( p.3[ ) when Tg = r 
and qgyg ^ 

Xg — {Xg, Xg) = [Xg, Qg H ~S:g) = 

P°s 

= {xl- ^ + (T, + ^^)kg) = xf-)'^ + e^(«(p.))a", 
= (x- Yg) = 

i~o I ^sys ^ Ps -1 . /rr , ^s ' Zg ,r> , 

= [Xg-, -[zg + ' + [Tg + A;, = 

eg eg[l + u°{pg)\ eg 

~^l I ^At {Ssys X Ps^'^ 



''^^'''^e,[l + <p,)] 
xf-)'^ + e^:(w(p,))a^. 



R's = (^:; Rs) 



^^'6.^^^ e.n°(p.)[l + <p,)]J + ^'^+ eg 

~^l _ i^sys XPsY _ 

Xg I/. — 



egu°{pg)[l + u°{pg)] 
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= X 



- e iuiv ))b(^-)'' - - 'rMPs))K{Ps)ST + ST] 

~ ^sQsvs + 



es[l + u"(ps)] 



~r . SlygU^{ps) _ r r rt 



_ _ Sl^gU^Ps) _ , , [1 - ^ 



U-{ps)[l+U-{ps)] ' 

^ -)'^(t) = /or g-;,,R^O, (3.7) 

namely in the gauge qgys ~ the external Fokker-Pryce non-canonical center of inertia 
coincides with the origin a;^^°'''^^(r) carrying the "internal" center of mass (coinciding with 
the "internal" MoUer center of energy and with the "internal" Fokker-Pryce center of inertia) 
and also being the Pirani centroid and the Tulczyjew centroid. 

Therefore, if we would find the center-of-mass canonical basis, then, in the gauge qgyg ~ 
and Tg pa r, the perfect fiuid configurations would have the four-momentum density peaked 
on the worldline x'f^^'^^{Ts); the canonical variables al^iir, a), Ureiiir, a) would characterize 
the relative motions with respect to the "monopole" configuration describing the center of 
mass of the fiuid configuration. The "monopole" configurations would be identified by the 
vanishing of the relative variables. 

Remember that the canonical center of mass lies in between the MoUer center of energy 
and the Fokker-Pryce center of inertia and that the non-covariance region around the Fokker- 
Pryce 4- vector extends to a worldtube with radius (the MoUer radius) \Ssys\/P'^- 
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IV. DIXON'S MULTIPOLES IN MINKOWSKI SPACETIME. 



Let us now look at other properties of a perfect fluid configuration on the Wigner hy- 
perplanes, always using the dust as an explicit example. To identify which kind of collective 
variables describe the center of mass of a fluid configuration let us consider it as a relativistic 



extended body and let us study its energy-momentum tensor and its Dixon multipoles ||20 
in Minkowski spacetime. 

The Euler-Lagrange equations from the action ( |1.10| ) for the dust] are 



(—-a 



r,(T 



0, 



(4.1) 



where we introduced the energy-momentum tensor [with a different sign with respect to the 
standard convention to conform with Ref. illll 



2 6S 



- P9 



, ( Ab 



B 



r, 0- 



A tB 



(4.2) 



When d^[y/gzg] = 0, as it happens on the Wigner hyperplanes in the gauge T, — r ^ 0, 
A(r) = 0, we get the conservation of the energy-momentum tensor T^^, i.e. dj^T^^ = 0. 
Otherwise, there is compensation coming from the dynamics of the surface. 

As shown in Eg. (|A9[) the conserved, manifestly Lorentz covariant energy-momentum 



tensor of the perfect fluid with equation of state p = p{n, s) [so that p = n-^\s — p)] is 



T^%x)[&] 



_.,v^ ^^^^^^^^ 

-e{p + p)UW +p^g^"']ix) 



dp 



dust 



-e/i 



(x). 



z"",.]^ = z^.r + '4f. (4.3) 



"A 

Therefore, in S^-adapted coordinates on each E^- we get 
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T^^(r, a)[a] = z^^{t, a)z^{T, a)T^'^{x = z{t, a))[a] - 



On Wigner hyperplanes, where Eqs.( p.38D hold and where we have 
z^{T,a) =x^(r) + ej:(w(p,))a", 



r=l 



get [A = A] 



n 



N 



is ■ u{ps) 



N 



dust 



pr 



T^''[^'s{r) + ei{u{ps))a-][a] 



= [S\i's{r) + S:,e>^HPs))WBi:{T) + 6%e:{u{ps))]T^''{r, 
= x^(r)xr(r)T-(r,a) + e^f(«(p,))6:(«(p,))T-(r,a) + 
+ [i^iT)e:{u{ps))+i:{T)e^{u{ps)W-{T,a), 



+ n 

dust 



[is ■ u{ps)f 
dp e 



+ 



(J 



T\2 



dn ^ [xs ■ u{ps)Y [is ■ u{ps)f[{rY - 6^,Y-Y-] 



— e 



[is ■ u{psy[' 



r, cr 



p^ + 5„,(T-i)«(T-i)-in,n, 

S'^is^ef^^Hps)) 



n 



+ 



dust 



[is ■ u{ps)f 
dn [xs ■ uips)]"^ 

j\r 

[is ■ u{ps)Y[{J-Y - 6uuY-Y 

r 



r, cr 



[is ■ u{ps)\'' 



p^ + SUT-'YKT-'Y^'n.U, (r,cr) 



,isf.e'i{uiPs))isue'^iu{ps)) . 



[is ■ u{ps)f 
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+ 



+ 



dust 

— —e 



Xs ■ uipsWixs ■ uipsWiiJ-y - 6uvY-Y- 

V/^2 + <5„„(T-i)-(r-i)-in,nJ (r, a) . 



[Xs ■ u{ps)Y J'' 



(4.6) 



Since we have 



x^,{r) = -A'^(t) = e[u^{p,X{p,) - e'^{u{p,))e';{u{p,))]xUr) 



— € 



u''{Ps)X{T)+e'^{u{ps))Xr{T) 



xl{T)^X'iT)-\{T)>0, 

U^^r) = ^^(^^ = ^-A(t)«'^(p.) + A.(T)6(f(«(p,)) 



x2(r) 



A2(r) - A2(r) 



(4.7) 



the timehke worldhne described by the origin of the Wigner hyperplane is arbitrary 
(i.e. gauge dependent): x'^{t) may be any covariant non-canonical centroid. As al- 
ready said the real "external" center of mass is the canonical non-covariant x'^{Ts) = 



x^{Ts)- 



PsuS^f'+esiS^Z+S: 



OV PsvVs 



: it describes a decoupled point particle observer. 



In the gauge Tg-r ^0, Xsys = Qsys ~ 0, implying A(r) = -1, A(r) = [grr = e, N = 1, 
— grr — 0], we get xi^iTs) — u^{ps)- Therefore, in this gauge, we have the centroid 



x':{Ts) = x^^^y^^^iTs) = x>:{0) + Tsuf^ips), 

— * 

which carries the fluid "internal" collective variable Xgyg — Qsys ~ 0. 

In this gauge we get the following form of the energy-momentum tensor [Y^ — J^] 

r'^1a;f-^^(T,) + e(^(«(p,))a«][a] = w^(p,)w^(j9,)T-(T„ <?) + 

+ [u>^ips)e':.{uiPs))+u''ipsKiu{ps)W^{Ts,a) + 
+ e'^{u{ps))e:{u{psj)r-'{Ts,a), 



(4.8) 



dp 



T\2 



+ n—\s [e - 



in 



dust 



—e 



T'^iTs, a) 



JV/i2 + 5,,(T-i)»(T-i)-m,n,J (T„ a). 



dp 
— n-^\s 



(Ts^a) 



dust 



-e [jV/i2 + <5„„(T-i)»(T-i)-m,n,] (T„ a), 

+ 



)](Ts,a) 
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with total 4 — momentum 
P^[a] = J d'aT'^''[x^{T,) + e^MPs))cr''][a]uM = 



and total mass 



M[a] = P^[a]uM = -P" = -Msys. (4.9) 

The stress tensor of the perfect fluid configuration on the Wigner hyperplanes is T^'^{Ts, a). 

We can rewrite the energy-momentum tensor in such a way that it acquires a form 
reminiscent of the energy-momentum tensor of an ideal relativistic fluid as seen from a local 
observer at rest (see the Eckart decomposition in Appendix B): 

T^''[a]=[p[a,U] u^{p^)u''{p^) + 

+ P[«,n] [r^^^- -u^{p^)u''{ps)] + 
+ u^{ps)q''[a, n] + n'^(p.)g^[a, U] + 
+ TZ[a,U] 6^f(n(p,))er(n(p,))](T„a), 

P[a,n] = ^^T- 

"J u 

q^[a,U] = e^Auip.))T^\ 
1 

— ( 
3 



T-[a, n] = - Y: 5™Tr[a, U] = 0, (4.10) 



where 

i) the constant normal u^{ps) to the Wigner hyperplanes replaces the hydrodynamic 
velocity field of the fluid; 

ii) p[a,U]{Ts,a) is the energy density; 

iii) V[a,Il]{Ts,a) is the analogue of the pressure (sum of the thermodynamical pressure 
and of the non-equilibrium bulk stress or viscous pressure); 

iv) q^[a, H](Ts, a) is the analogue of the heat flow; 

v) T^^[a,I\\{Ts,a) is the shear (or anisotropic) stress tensor. 



We can now study the manifestly Lorentz covariant Dixon multipoles pO[ for the per- 
fect fluid configurationon the Wigner hyperplanes in the gauge A(r) = —1, A(r) = [so 
that x^(T,) = u^'(p,), i^(T,) = 0, x^(T,) = xf^y^^^'{T,) = a;^(0) + u''{ps)Ts] with re- 
spect to the origin an arbitrary timelike worldline w^iTg) = z^^{Ts^fj{Ts)) = x^f^i^'^^^iTs) + 
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e!fiu{p,W{Ts). Sincewehavez'^(r„a) = x'f^y^^^'{T,)+e>f{uips))a^ = w^{T,)+e>^{uips))[a-- 
V'^i^s)] w'^{Ts) + Sz'^{Ts, a) [for f]{Ts) = we get the multipoles with respect to the origin 
of coordinates], we obtain [ means symmetrization, while [//i../^„] means antisym- 

metrization; t^i-^"'^"(T„ rf = 0) = t^i-^"'^"(T,)] 

+ [u^{PsK.{u{ps)) + «'^(p.)elf (tz(p.))]/;-^"^-(T„ rf); , 
/^-'■"^^(r„77) = I dVla'-^ -7y^Hr,)]...[a'-" -ry'-"(r,)]r^^(T„a)[a], 

For?/ - n = {monopole) I^^{T,) = P\ It{Ts) = P\ 
4'-''-%{Ts) = J c^V5xf (a)r%[xf-)'^(r,) + 6^(u(p,))a"][a] = 

The Wigner covariant multipoles I^'-''"^^{Ts), I^'-'^''{Ts), I^'-'^^''^ (T,) are the mass. 
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stress and momentum multipoles respectively. 

The quantities I^^'^^^a{Ts) and i^^'^{Ts) are the traceless quadrupole moment and the 
inertia tensor defined by Thorne in Ref. ||21|| . 

The quantities /^^'"^"^"^(T^) and i^^'^(Ts) are Dixon's definitions of quadrupole moment 
and of tensor of inertia respectively. 

Moreover, Dixon's definition of "center of mass" of an extended object is tj^^{Ts) = 
or I'^'^iTs) = —P^rlyg = 0: therefore the quantity fsys defined in the previous equation is 
a non-canonical [{r^^^, rfj^^} = S^yg] candidate for the "internal" center of mass of the field 

configuration: its vanishing is a gauge fixing for P ^ and implies ^^(T,) = x^^''y''''^{Ts) = 
x^(0) + u^{j)s)Ts- As we have seen in the previous Section Tgyg is the "internal" M0ller 
3-center of energy and we have Tsys ~ qsys ~ Vsys- 

When Fjf^iTs) = 0, the equations = ^^^^^ = -P^^ = ^ - implies the correct 

momentum- velocity relation = ^^^^ ^ 0. 
Then there are the related Dixon multipoles 

n = ^ p^(T,) = P^[a\ = -ee'^{u{ps))P^ ^ -ep^, 

pr''"'{Ts)uM = t^^-^"(T,) = e^f;(«(p.))...etf:(n(p.))j;— -(T,). (4.12) 
The spin dipole is defined as 

5r(T.)H = 2p^r\T,) = 2e\^{u{ps))ei{u{ps))r/^{%) = 

= = 

= e'^{u{ps))e':{u{ps))S:^y, + [e^{u{ps))u''{ps) - e';{u{ps))u^{ps)]Sl^s, 

w^(p.)5rm)N = -eU<Ps))S:^s = -iUTs) = P^e:.{u{Ps)ysys, (4.13) 

with Ufj,{ps)Sjy (Ts)[a] = when ii^{Ts) = and this condition can be taken as a def- 
inition of center of mass equivalent to Dixon's one. When this condition holds, the 
barycentric spin dipole is Sj!^{Ts)[a] = 2eli^{u{ps))e'^^{u{ps))I^^'^{Ts), so that I^^^'^{Ts) = 



As shown in Ref. |20|, if the fiuid configuration has a compact support W on the Wigner 
hyperplanes T^wt and if f{x) is a complex- valued scalar function on Minkowski spacetime 
with compact support [so that its Fourier transform f{k) = J d!^xf{x)e^^'^ is a slowly increas- 
ing entire analytic function on Minkowski spacetime {\{x° + iy''Y° ...{x^ + iy^Y^ f{x^ + iy^)\ < 
C<jo...g3e"°'^°'''''"'''"^'^ ', > 0, positive integers for every /x and Cg^,,,g^ > 0), whose inverse 
is f{x) = f (|^/(fc)e~*'^''^], we have [we consider ff = with 5z^^ = 5x'^] 
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< T>"'J> = J d^xT^''{x)f{x) = 
dT, f d^a 



oo (_j\n 



n=0 

dT. 



nl 



(27r) „=o 
and, but only for f{x) analytic in W we get 

r °° 1 



(4.14) 



°° 9" fix 

n=0 



4 



T^^(x)[a] = 5: 



-1)" 



/t — u 

6:;(n(p.))6^(«b.))T^^(T„a)[«] 



n=0 



J7=0- 



(4.15) 



For a non analytic /(x) we have 



TV ^ 



(9" f fx) 



+ dT, 



n=N+l 



nl 



i^fii---i^tj.„ i^T 



(Ts), (4.16) 



and, as shown in Ref. |2^, from the knowledge of the moments tj^"''^"'^{Ts) for all n > we 
can get T^^{x) and, thus, all the moments with n < N . 

In Appendix D other types of Dixon's multipoles are analyzed. From this study it turns 
out that the multipolar expansion (|4.15|) may be rearranged with the help of the Hamilton 
equations implying d^j^T^^ = 0, so that for analytic fluid configurations from Eq.( |D5|) we get 

T^''{x)[a] ^ -eM(^(p,)e^(M(p,)) j dTj\x - XsiT,)) + 
1 d 



+ 



2dxP 



n=2 



J dTs 5\x - xs{%)) S!^^''{T,)[a] u^\ps) + 
-1)" 9" 



n\ dx^^^ ...dx'^" 



rfT,54(x-x,(T,)) J^-^"^^(r,), 
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n 

+ 4'-'-''''\T„ff)]\, (4.17) 

where forn > 2 and r/ = X!^^-^-^\T,) = ^i^J(!^^■■^^n-lM^^n)u^J.^^^ ^^^^ quantities 
jjj.L.tj.ntJ.i'pcr ^rp^^ being the Dixon 2^+"-pole inertial moment tensors given in Eqs. (P7|) [the 
quadrupole and related inertia tensor are proportional to I^^^^'^'^iTg)]. 

The equations d^T^^ = imply the Papapetrou-Dixon-Souriau equations for the 'pole- 
dipole' system P^{Ts) and S^''{Ts)[a] [see Eqs.(|nT]) and (|D|); here ff = 0] 

^^£^M^2P),^{Ts) u''\ps) = -2eP^e^nu{Ps)) W'Kps) ^ 0. (4.18) 

The Cartesian Dixon's multipoles could be re-expressed in terms of either spherical or 
STF (symmetric tracefree) multipoles |^T| [both kinds of tensors are associated with the 
irreducible representations of the rotation group: one such multipole of order / has exactly 
2/ + 1 independent components] . 
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V. ISENTROPIC AND NON-ISENTROPIC FLUIDS. 



Let us now consider isentropic (s = const.) perfect fluids. For them we have from 
Eqs.( |1.9| ), ( p. .101 ) and (|2.1| ) [in general /i is not the chemical potential but only a parameter] 



n 



L = -^N^fi^). (5.1) 

Some possible equations of state for such fluids are (see also Appendix A): 
1) p = 0, dust: this implies 

p{n) = pn = P~^-' 

2) p = kp{n) = n^^^ — p{n) {k ^ —1 because otherwise p = const., p = 0, f = —Ts). 
For = I one has the photon gas. The previous differential equation for p{n) implies 

p(n) = (an)'=+i = /in^+^ = p{—)''^\ {p = aJ^^^), i.e. 

(5,, 

V7 oX ^ ^ 



[for A; ^ we recover 1)]. 

More in general one can have k = k{s): this is a non-isentropic perfect fluid with p 
p{n, s). 

S) p = kp^{n) = n^jj^ — p{n) (7 7^ 1) [^. It is an isentropic polytropic perfect fluid 



dn 

(7 = 1 + -). The differential equation for p{n) implies [a is an integration constant; the 
chemical potential is = 



an an 

n\Ti} = = Z C 

[1 - A;(an)^-i]^ [1 - A;(an)^]" ' 



X X 



7 [l_A;(a^)7-i]^ [1 _ A:(a;|)^]n 



= —[1 - Ma— )^-r^ = —[1 - Ka—)^-^-^'\ (5.4) 
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Instead in Ref. p3| , |2^ a polytropic perfect fluid is deflned by the equation of state [see 
the last of Eqs. (p^ ); m is a mass] 

k(s) 

p(n,ns) = mn -\ (mnV, (5-5) 

7-1 

and has pressure p = k{s){mny' = (7 — l)(p — mn) and chemical potential or specific 
enthalpy fi = mc^ + mA;(s);^(mr2)'^~^ of Eq.(^ 

4) p = p{p), barotropic perfect fluid. In the isentropic case one gets p = p{n) by solving 
p(p(n)) = p{n). 

5) Relativistic ideal (Boltzmann) gas p| (this is a non-isentropic case): 

p = nksT and p = mc'^nT(P) — p, p = ^ ^ = mc^r(/3), (5.6) 

with P = -^T, T(P) = {Ki are modifled Bessel functions). One gets the equation of 

state p = p{n, s) by solving the differential equation 

^1 -mc'T(^^^^) (5 7) 

5a) Ultrarelativistic case « 1 (mc^ « ksT): since we have r(/3) ^ 1 + f + 0(/53)^ 
we get p = SnksT + + 0{kBTl3^), namely p ~ |p and p ^ p{n) = pn^^^. 

5b) Non-relativistic case (3 » 1 [ksT « mc^): since we have r(/5) ~ 1 + ^ +0(/3~^), 

we get p ~ mc^n+|p, so that we have to solve the differential equation 3n|^|s— 5p+2mc^n ^ 
0. Its solution is p{n,s) ~ mc^n + k{s)n^^^. To flnd k{s) let us use the deflnition of 
temperature: p = nksT = A;_b|j|„ = ksn^^^^^^ = ''^f^ls ~ P = ^n^^^k{s). This leads to the 

equation dlnk{s) = whose solution is k{s) = he^^B = e with h = e '■^^b = const.. 
Therefore, in this case we get [it is a polytropic like in Eq.( p.5D with 7 = 5/3 and k{s) = 

2(s-So) 

, , 2 ^ 2(s-so) I dp, 2 5 

p{n, s) ^ mc n + n^e ^''b and T=— — \n^—. — n^e ^''b . (5.8) 



n ds 3k 



B 



The action for these fluids is [s = s{a^)] 

S = J dTd^aL{a\T, a), z^{t, ^)) = - J drd^aN^piji, s), (5.9) 
and we have as in Section II 



3 



k 



1=1 



h 3i = cychc, 



J| = JN^iJ^y - ^gf,[r + mr][J^ + N-"r] = NX, 
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Y 



N 



dX Y^Tfi 



ddrai NX ' 

dX _ ^gfsY'Y' _ {rf - X'^ 



{i,j,k cyclic), 



dX 



-J 



NX 

'QusY' 



NX 



NX 



(5.10) 



In the cases 1), 2) and 3) [in case 3) we rename fj, the constant a] the canonical momenta 
can be written in the form 



n,(r, a) 



dL{T, a) 



ddTa}{T, a) 



rdf{x) 


\rfrp 
1 J. 


dx 


X \ 



dx 



^ i rdfa' N'r-NY\ , 
a-ra — — = Of-a 



= N'dfa' + Ndfa\T-y%^[^^ 



dL{T, a) 



I 



III. „ - . t^s/i y 



"ii 

-X 



liJ^' dx '^=v^^ ' 
dL 



" dN " dN^ 



dx 



+ 



X dx '^=^ 



X 

-e[r{T-yil,Zr,]{T,a) 



with 



Vt 7 



X 

V7 



dx 



2L = f{n) + 



^ - 5/(n) 



n 



dn 



(5.11) 



To get the Hamiltonian expression of the constraints H'^{r, a) ft! 0, we have to find the 
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solution X of the equation + ^gf^Y^{X)Y''{X) = {J^f with Y^{X) given by the second 
hne of Eq. (|5.11|) . This equation may be written in the following forms 



n 



X' 



ox 



\-2 



1 + 
1 + 



df{x) 



^2 _|_ ^2 V Qfi 



\-2 



5", or 
- 1 



^2 = ^g,,{T-^Ym,{T-ym 
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2/ Tr\2 



52 



52 



or 



^{ij? + A^y 



A2 I JT\2 



7 

^2 52 



52 



/x2 + A2'7(/i2 + A2)- 



/i2 + A2'7(/i2 + ^2) 



ML + A^^z, 



(5.12) 



Therefore, all the dependence on the metric and on the Lagrangian coordinates and 
their momenta is concentrated in the 3 functions ^7, A? j ^ = ^gfs{T~^)"Ili{T~^)'^^Ilj/ fi^, 

S2/^2^ = (J-)2/^. 

Let us consider various cases. 

1) p = 0, dust. As in Section II the equation for X and the constraints are 

X2 [/x2 + A2]=52, 



X 



B 



y/J^^TA^ .//i2 + 3^._(T-i)"n,(T-i)^m, 



>A^^o\r\[l-^ + 0{A% 



\.r\{T-^Ym. 



fi^ + ^grs{T-'ym,{T-^ym, 



2) p = kp, k ^ —1. The equation for X is 

^2 



X2 [^2 



(A: + l)2(^)2/cJ 



5^ or 



1 



fA)2 Ml { 

V ^p2 + ^2l(^ + l)2(^)2fc 



1)] 



fi2 



7(p2 + A2)^ 



(5.13) 
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Let us define Z as the deviation of X from dust (for — > [S^a* = 0]: we have 



r 



—JSl — Z. Then we get the following equation for Z 



P' 



+ 



ZV + PlZ-'"]^!, 



1, or 



a 



Pi 



P" 



+ 

A2(//2 + A2)fc-l^fc 

(A; + 1)2^2^ 



We may consider the following subcases: 
2a) k — m ^ —1, with the equation 



151 



Zi 



[a2 + /3^Zr-] = l, or + a^Zi - 1 = 0. 

i) p = p (/c = m = 1), with the equation 



VS2 



./4^ - >12 ^ 
X = ^^^^^ = — W4/x^ 



2^^ + 3^,,(r-i)-n,(r-i)^in,-, /or ^2<4__ 

7 



ii) p = 2p (A; = m = 2), with the equation 



a^Zl - Zi + /32^ = 0, Zi 



1 



2q;2' 



[l±Vl-4a2/?2], 



iii) p — —2p {k — m — —2), with the equation 

(3l^Zl + a^Zi -1 = 0, 
2b) k——, with the equation 
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(5.20) 



i) P = = |) ^ = 2), with the equation 



1 



ii) p — —\p {k — — |, m — —2), with the equation 

iii) p = photon gas (A; = |, m = 3), /5^/3 = ieg2/'3^^l!,.^^2)2/3 , with the equation 

a^Zl + l3l/^Zl-l = Q, or Zl+pZ^,+r = 0, 



(5.21) 



(5.22) 



Z2 
X 



with p 

1 
3 
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/5?/3 1 

^ > 0, r = ^ < 0, 



Yo P = -^2 ^ 

3^ 3q;2 



1/3. 



_ 3^4V+i!W^^3/2 

'^^ ^ 16^2/3 > ' 



P 



1/3 



27.4i(//- + A- 



2\2/3^ 2/3 



3q;4 



28^4/3 



(2^-2r)-^( 



— <0, 

+ ^2 / 27^67 



4 ^ 27 



1 



1 • 27tt^ 



(27-4 



/^l/3^ 



(/i2 + A2)2^^_ 27^67 



1)' 



4/x'i 



2iV^5^ 



F2 



X 



3^271/3 



+ 



1 / 



16(J-)2/3 2V3V 2^^ii^{J-f \ 2iV'(^") 
27^67 



27^67 -^1/3 



2V3 



2iV'(^")' \ 2iV'(^^)' 

9^2^1/3 



27^67 xl/3n 



3/2 



32(J^)2/3 



+ 0(A^)], 
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4 



2^. 



(5.23) 



In the case of the photon gas we get a closed analytical form for the constraints. 
3) p = kp^ , 7 = 1 + -, 7 7^ 1 7^ 0), with the equation 



fx' 



7 
X 



2(n+l) 



or 



or 



V7^ 



1 + 



/i2 + A2 



B"' 



7(/i2 + ^2)' 



^ a/7 



= pV7 



(ZD! 

7 



(5.24) 



Let us define Z as the deviation of X from dust (for Q [drOt^ = 0]- "we have Z — > 1, 

X = J'^^ Z). Then we get the following equation for Z 



P^ 



/i2 + 



+ 



(A; + 1)25 



2fc 



fi\B\Z 



IN 2(n+l) 



y7v/p^ + ^' 



(5.25) 



In conclusion, only in the cases of the dust and of the photon gas we get the closed 
analytic form of the constraints ( i.e. of the density of invariant mass 7Vl(r, a), because for 
the momentum density we have M."^ = J'^ {T~^)"Ili independently from the type of perfect 
fluid). 

In all the other cases we have only an implicit form for them depending on the solution 
X of Eq. ( |5.12|) and numerical methods should be used. 
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VI. COUPLING TO ADM METRIC AND TETRAD GRAVITY. 



Let us now assume to have a globally hyperbolic, asymptotically flat at spatial infinity 
spacetime with the spacelike leaves of the foliations associated with its 3+1 splittings, 
diffeomorphic to |2|^|lTl. 

In (TT--adapted coordinates a"^ = (o""^ = r; a) corresponding to a holonomic basis [da^, 
Oa = d/da"^] for tensor fields we have [A^ and A^'' are the lapse and shift functions; ^grs 
is the 3-metric of S^; /^(r, a) is the unit normal vector fields to S^] 

= 9AB = {Vr = e(A^^ - VsA^'iV'); Vr = -e^g-rsN'; = -e Vs} = 

^AB = ^Qab - eUh- (6.1) 

A set of ST--adapted tetrad and cotetrad fields is (a) = (1), (2), (3); ^e^"^-, and ^e^^^ = ^e(^a)r 
are triad and cotriad fields on T,^] 

4 T^A _ lA _ f ^ . \ 

(s)^(o)-« ~^iV'~lV^' 
f^^E^°^ =Ia = {N;0), 

f^)E^^^ = (iV(») = iV 3e("); =^e('^)). (6.2) 
In these coordinates the energy-momentum tensor of the perfect fiuid is 

T^^ = -e(p + p)f/^f/^+pV''- (6.3) 
If we use the notation F = e/^f/^ = eNU"^ , we get 

Tab = ''gAc'^gBcT^''' = [dJc + AacWd + A^^jT^^ = 
= EIaIb + JaIb + IaJb + Sab, 

E = T^''lAlB = -4ip + P)^'-pl 
J A = eAAcT^'^'lB = -e{p + p)TAABU'', 
Sab = Aac^bdT'''' = -t[{p + p)AacU'' AbdU'' - pAab]- (6.4) 

The same decomposition can be referred to a non-holonomic basis [9^ = {6^ = Ndr; 9^ = 
da' + N^'dr}, Xa = {Xi = N-\dr - A^'^,); dr}; A = (/; r)] in which we have [f = ehU^ = 



eU^ = yl — ^grsU^'U'^, so that W may be interpreted as the generalized boost velocity of 
with respect to P [11,11; ^ 



9Ab = { 9ii = e; 9ir = 0; Qrs = 9rs = -e 9rs} 
= elAls + ^AB, 
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All = A;^ = 0, A^, = -e^Qrs: ^AbU^ = (0; -e^grs): 

Tab = -^(p + p)UaUb + p'^gAB = 

— EI^Ib + Ja^B + ^aJb + ^ABi 



E^f^%lB = -e[{p + p)r'-p], 
-jA = eA^c2"^^^5 = Qi = 0; = -e(p + p)f 
^Ac^bdT'^'^ — 



(6.5) 



and jV are the energy and momentum densities determined by the Eulerian observers 
on while S^s is called the "spatial stress tensor". 

The non-holonomic basis is used to get the 3+1 decomposition (projection normal and 
parallel to S^) of Einstein's equations with matter "^G^^ = ■^qT'^^ [when one does not has 
an action principle for matter, one cannot use the Hamiltonian ADM formalism]: in this 
way one gets four restrictions on the Cauchy data [ ^G" = ^g-^'' '^Qir ^ j£_f'ir.^ ^^^ey 
become the secondary first class superhamiltonian and supermomentum constraints in the 
ADM theory; k = c^/8nG] 



o 2 - 
[r, d) = -E{T, d), 



(6.6) 



and the spatial Einstein's equations ^G^^ = ■^qT^'^ ■ By introducing the extrinsic curvature 
^Krsi this last equations are written in a first order form [it corresponds to the Hamilton 
equations of the ADM theory for ^Qrs and — ■^^/l^K''^ — ^g'" '^K); "|" denotes the 
covariant 3-derivative] 



dr ^Krs(r, a) = {N[^Rrs + ^K^Krs - 2 ^K^s\ - 

- N\s\r + iV"|, ""Kur + N^\r ^K^s + ^Krs\u) (t, d) 



k 



{Srs + ^^9rs{E-S\)\{T, d) . 



(6.7) 



The matter equations T^^.^^, = become a generalized continuity equation [entropy con- 
servation when the particle number conservation law is added to the system] JaT^^-b ^ 
and generahzed Euler equations A^bT^'^-c ^ [^x is the Lie derivative with respect to the 
vector field X] 



drE + Nf\r = N{S'' -"Krs + E K) - 2fN\r + Cf^E 



dijr + NS^'\s] (r, d) ^ \n{2 'K^% + f 'K) - S^'N\, - EN^^ + C^f] (r, d). 



(6.8) 
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The equation for Srs would follow from an equation of state or dynamical equation of 
the sources [for perfect fluids it is the particle number conservation]. 

This formulation is the starting point of many approaches to the post-Newtonian ap- 



proximation (see for instance Refs. p3| , |30[| ) and to numerical gravity (see for instance Refs. 

mm)- 

Instead with the action principle for the perfect fluid described with Lagrangian coordi- 
nates (containing the information on the equation of state and on the particle number and 



entropy conservations) coupled to the action for tetrad gravity of Ref. ||2^ (it is the ADM 



action of metric gravity re-expressed in terms of a new parametrization of tetrad fields) we 
get 



S = -ek I drd^a {N e(a){6){c) ^e[„) ^e^^) '^^rs{c) + 



3^ 



e 

2iV' 

- / dTd'a{N^p{^, .)}(r, a). (6.9) 



The superhamiltonian and supermomentum constraints of Ref. p5[ are modified in the 
following way by the presence of the perfect fiuid 

nr^Qr = Qor + Mr^O, (6.10) 

with Mr = J^iT-^ym, = -drami and with M given by Eq.(|2J^) for the dust and by 
Eq.( [5.23D for the photon gas. 



In the case of dust the explicit Hamiltonian form of the energy and momentum densities 

is 



Mir, a) = rir,a)Jf,^+[^grsiT-^ia)Y^iT-^{a)ym,U,]iT,a) 



-.et(..«^(.,.)),/,^ + 3,..^^n.n. 



Mrir,a) = 'grsriT,a)[{T-\a)rmr,a) = -9,a^(r, a)n,(r, a). (6.11) 

In any case, all the dependence of Ai and Ai^ on the metric and on the La- 
grangian coordinates and their momenta is concentrated in the 3 functions ^^7, A^/fj!^ = 



The study of the canonical reduction to the 3-orthogonal gauges |26|,0 will be done in 
a future paper. 
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VII. NON-DISSIPATIVE ELASTIC MATERIALS. 



With the same formahsm we may describe relativistic continuum mechanics [any rela- 
tivistic material (non-homogeneous, pre-stressed,...) in the non-dissipative regime] and in 
particular a relativistic elastic continuum P [see also Refs. p3| , |3^ and their bibliography] 
in the rest-frame instant form of dynamics. 

Now the scalar fields a*(z(r, a)) = a*(r, a) describe the idealized "molecules" of the 
material in an abstract 3-dimensional manifold called the "material space", while J"^(r, a) = 
[N-^nU^]{T, a) is the matter number current with future-oriented timelike 4- velocity vector 
field U^{t, (t); n is a scalar field describing the local rest-frame matter number density. The 
quantity 9^a*(r, a) = z^(r, cr)(9^a*(z) is called the "relativistic deformation gradient" in 
ST--adapted coordinates. 

The material space inherits a Riemannian (symmetric and positive definite) 3-metric 
from the spacetime 

Qio = ^gi^-d^ctd,cP = ^g^^dAa'd^a^. (7.1) 

Its inverse Gij carries the information about the actual distances of adjacent molecules 
in the local rest frame. 

For an ideal fluid the 3-form rj of Eq( |1.3| ) gives the volume element in the material 
space, which is sufficient to describe the mechanical properties of an ideal fluid. Since 
we have that n = 4^ is a scalar, we can evaluate n in the local rest frame at z where 

U^{z) = ( , ^ _ ; ) and do<y^\z = [in the local adapted non-holonomic basis we have 
V^o/V^ = VW^- We get at z 



n = = 77123 det dka' ^ = 17123 det dka' \fdet\^g^\ = 77123 V^et G'^ . (7.2) 



Therefore, we have n = rjusVdet G^K 

Moreover, the material space of elastic materials, which have not only volume rigidity but 
also shape rigidity, is equipped with a Riemannian (symmetric and definite positive) 3-metric 
7|*^''(a*), the "material metric", which is frozen in the material and it is not a dynamical 
object of the theory. It describes the "would be" local rest-frame space distance between 
neighbouring "molecules" , measured in the locally relaxed state of the material. To measure 
the components '~)\j ^\a'^) we have to relax the material at different points a*(r, a) sepa- 
rately, since global relaxation of the material may not be possible [the material space may 
not be isometric with any 3-dimensional subspace of M^, as in classical non-linear elas- 
tomechanics, when the material exhibits internal stresses frozen in it]. The components 
7^*^'' = 7f*^''(a*(r, a)) are given functions, which describe axiomatically the properties of 
the material (the theory is fully invariant with respect to reparametrizations of the material 
space) . 

For ideal fluids 7ij^'' = ^if this also holds for non-pre-stressed materials without "inter- 
nal" or "frozen" stresses. 
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Now the material space volume element t] has 

^i23(«*) = \J d^'^ ) SO that n = 77123 V det G*^ = det y/ det G^^ , (7.3) 

and it cannot be put = 1 like for perfect fluids. 

The pull-back of the material metric 7^/^'* to M'^ is 

7?fi = l^^d^d^cP satisfying ^fj^U^ = 0. (7.4) 

The next step is to define a measure of the difference between the induced 3-metric 
Gij{da^) and the constitutive metric '~)\^\a^)^ to be taken as a measure of the deformation 
of the material and as a definition of a "relativistic strain tensor", locally vanishing when 
there is a local relaxation of the material. Some existing proposal for such a tensor in 
are 



33,34 



S% = \i"9AB-^U^U^-^f^) (7.5) 



which vanishes at relax and satisfies S'^^j^U^ = [but it must satisfy the involved matrix 
inequality 2 det S^^^ > det {^g - eUU)]. 

u) S^^^^^ = liK^^-5i), ^(^)^^f/^ = 0, (7.6) 

with K^j^ = ^9^^{l^cS ~ ^^cUb), K^bU^ = [the 4- velocity field is an eigenvector of 
the i^-matrix]. 

in) S^^^ = -hnK, (7.7) 

with the same i^-matrix as in ii). 

However a simpler proposal [^ is to define a "relativistic strain tensor" in the material 
space 

S^' = llk^G''^ (7.8) 

with locally Si^ = 6l when there is local relaxation of the material (in this case physical 
spacelike distances between material points near a point z'^{T,(f) agree with their material 

distances). Since n = 77123 Vrfet G^^ = \J det 'yl^^^ V det G^^ , we have n = V det Si^ for the local 
rest-frame matter number density. 

We can now define the local rest-frame energy per unit volume of the material 
rz(r, ?)e(r, a), where e denotes the molar local rest-frame energy (moles = number of parti- 
cles) 

e(r, a) = m + u/(r, a). (7.9) 

Here m is the molar local rest mass, uj is the amount of internal energy (per mole of 
the material) of the elastic deformations, accumulated in an infinitesimal portion during the 
deformation from the locally relaxed state to the actual state of strain. 
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For isotropic media ui may depend on the deformation only via the invariants of the 
strain tensor. 

Let us notice that for an anisotropic material (like a crystal) the energy ui may depend 
upon the orientation of the deformation with respect to a specific axis, reflecting the mi- 
croscopic composition of the material: this information may be encoded in a vector field 
i?*(r, a) in the material space and one may assume ui = ui{G~j^E^E^). 

The function e = e[a\ ■■■] describes the dependence of the energy of the material 

upon its state of strain and plays the role of an "equation of state" or "constitutive equation" 
of the material. 

In the weak strain approximation of an isotropic elastic continuum (Hooke approxima- 
tion) the function ui depends only on the linear {h = Si'') and quadratic {q = Si^S^i) 
invariants of the strain tensor and coincides with the standard formula of linear elasticity 
[V = ^ = ^^^g J specific volume], 

Ui = \{V)h'^ + 2fi(y)q + 0{cubicinvariants), (7-10) 

where A and fi are the Lame coefficients. 

The action principle for this description of relativistic materials is 

S['^g,a\da'] = J dTd^aL{T,a) = - J dTd^a{N^){T, a) n{T,a) e{T, a). (7.11) 



It is shown in Refs. [|6|,|33| that the canonical stress-energy- momentum tensor Tg = 
pfd^a'' — 5^L, where pf = —-g§^ is the relativistic Piola-Kirchhoff momentum density, 
and coincides with the symmetric energy-momentum tensor T^^ = —2^^^, which satisfies 

T^^ .j^ = due to the Euler-Lagrange equations. This energy-momentum tensor may be 
written in the following form 

TAB = N^n[eUAB + ZABl (7-12) 

where = Zij d^a^dga^ is the pull-back from the material space to of the "response 
tensor" of the material 

de 1 ■• 

Zij. = 2-—, sothat de{G) = -ZijdG'^. (7.13) 
oG^-^ 2 

The part T^^j = nZ^^, r^^ f/^ = 0, may be called the relativistic "stress or Cauchy" 
tensor and contains the "stress-strain relation" through the dependence of Zij on Si^ implied 
by the consitutive equation e = e[a, 7*-*^'', G, ..] of the material. 

For an isotropic elastic material we get 



% = V 



pG-.^ + 57if + CG J , (7.14) 



where p = B = ^^ff, C = and we get 



de{V, h, q) = -pdV + ^VBdh + ^VCdq. (7. 15) 
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The response parameters describe the reaction of the material to the strain: p is the 
"isotropic stress" , while B and C give the anisotropic response as in non-relativistic elesticity 
[perfect fluids have e = eiV), B = C = 0, Zij = VpG^j^ and de{V) = —pdV is the Pascal 
law]. 



See Ref. [35| for a different description of relativistic Hooke law in linear elasticity: there 



is a 4-dimensional deformation tensor S^i, = I^V^^^+^V^^^) and the constitutive equations 
of the material are given in the form T^'^ = C^^^^'^"'^^ S^p- 

In Ref. the theory is also extended to the thermodynamics of isentropic flows (no 
heat conductivity). The function e is considered also as a function of entropy S = S{a) and 
de = ^ZijdG^^ is generalized to 

de = ^ZijdG'^ - SdT. (7.16) 

Then e is replaced with the Helmholtz free energy / = e — TS so to obtain 

df = ^Z,,dG'^ -SdT, (7.17) 

[for perfect fluids we get deiV, S) = —pdV + TdS, df(y,T) = —pdV — SdT]. This suggests 
to consider the temperature T as a strain and the entropy S as the corresponding stress and 
to introduce an extra scalar field a'^{T, a) so that T = const.U^dj^a^ . The potential a'^(r, a) 
has the microscopic interpretation as the retardation of the proper time of the molecules with 
respect to the physical time calculated over averaged spacetime trajectories of the ideahzed 
continuum material. 

In this case the action principle becomes S = — J drd^a Ny/jnf{G,T) {T,a) and one 
gets the conserved energy-momentum tensor T^^ = N^n[Q' + TS)U^U^ + Z^j^]. 

In all these cases one can develop the rest-frame instant form just in the same way as it 
was done in Section II and III for perfect fluids, even if it is not possible to obtain a closed 
form of the invariant mass. 
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VIII. CONCLUSIONS. 



In this paper we have studied the Hamiltonian description in Minkowski spacetime as- 
sociated with an action principle for perfect fluids with an equation of state of the form 
p = p{n,s) given in Ref. [||], in which the fluid is descrbed only in terms of Lagrangian 
coordinates. 

This action principle can be reformulated on arbitrary spacelike hypersurfaces embedded 
in Minkowski spacetime (covariant 3+1 splitting of Minkowski spacetime) along the lines of 



Refs. |TO| , pT| . At the Hamiltonian leve the canonical Hamiltonian vanishes and the theory 
is governed by four first class constraints ?^'^(r, a) ^ implying the independence of the 
description from the choice of the 3+1 splitting of Minkowski spacetime. 

These constraints can be obtained in closed form only for the 'dust' and for the 'photon 
gas'. For other types of perfect fluids one needs numerical calculations. After the inclusion 
of the coupling to the gravitational field one could begin to think to formulate Hamiltonian 
numerical gravity with only physical degrees of freedom and hyperbolic Hamilton equations 
for them [like the form ( p.50|) of the relativistic Euler equations for the dust]. 

After the canonical reduction to 3+1 splittings whose leaves are spacelike hyperplanes, we 
consider all the configurations of the perfect fluid whose conserved 4-momentum is timelike. 
For each of these configurations we can select the special foliation of Minkowski spacetime 
with spacelike hyperplanes orthogonal to the 4-momentum of the configuration, 

This gives rise to the "Wigner-covariant rest-frame instant form of dynamics" |TD|JTI| 
for the perfect fluids. After a discussion of the "external" and "internal" centers of mass 
and realizations of the Poincare algebra, rest-frame Dixon's Cartesian multipoles ||20| of the 
perfect fluid are studied. 

It is also shown that the formulation of non-dissipative elastic materials of Ref. [§, 
based on the use of Lagrangian coordinates, allows to get the rest-frame instant form for 
these materials too. 

Finally it is shown how to make the coupling to the gravitational field by giving the 
ADM action for the perfect fluid in tetrad gravity. Now it becomes possible to study the 
canonical reduction of tetrad gravity with the perfect fluids as matter along the lines of Refs. 

2M3. 
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APPENDIX A: RELATIVISTIC PERFECT FLUIDS. 



As in Ref. [|l|] let us consider a perfect fluid in a curved spacetime with unit 4- 
velocity vector field U^{z), Lagrangian coordinates a^{z), particle number density n{z), 
energy density p{z), entropy per particle s{z), pressure p{z), temperature T{z). Let J^{z) = 
\J^g{z)n{z)U'^{z) the densitized particle number flux vector field, so that we have n = 

\J £ '^Qiiv J'^ I V^- Other local thermodynamical variables are the chemical potential or 
specific enthalpy (the energy per particle required to inject a small amount of fluid into a 
fluid sample, keeping the sample volume and the entropy per particle s constant) 

P = -{p + p), (Al) 
n 

the physical free energy (the injection energy at a constant number density n and constant 
total entropy) 

a = - - Ts, (A2) 
n 

and the chemical free energy (the injection energy at constant volume and constant total 
entropy) 

f = -{p + P)-Ts = fi-Ts. (A3) 
n 

Since the local expression of the first law of thermodynamics is 

dp = fidn + nTds, or dp = ndfi — nTds, or d{na) = fdn — nsdT, (A4) 
an equation of state for a perfect fluid may be given in one of the following forms 

p = p{n,s), or p = p{p,s), or a = a{n,T). (A5) 
By definition, the stress-energy-momentum tensor for a perfect fluid is 

T^"" = -e pU^'U" + p^g^"^ - eWW") = -e{p + p)U>'U'' + p ^g^"", (A6) 
and its equations of motion are 

T^^, = 0, {nU^),, = -^Qp^J^ = 0. (A7) 

As shown in Ref. m an action functional for a perfect fluid depending upon J^{z), 
'^Ofiuiz), s{z) and a* (2;) requires the introduction of the following Lagrange multipliers to 
implement all the required properties: 

i) 0{z): it is a scalar field named 'thermasy'; it is interpreted as a potential for the fluid 
temperature T = In the Lagrangian it is interpreted as a Lagrange multiplier for 

implementing the "entropy exchange constraint" (sJ^)^^ = 0. 
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ii) ^p{z): it is a scalar field: it is interpreted as a potential for the chemical free energy 
/. In the Lagrangian it is interpreted as a Lagrange multipliers for the "particle number 



conservation constraint" J^^n = 0. 



iii) /3i{z): they are three scalar fields; in the Lagrangian they are interpreted as Lagrange 
multipliers for the "constraint" a^^f^J^ — that restricts the fiuid 4-velocity vector to be 
directed along the fiow lines = const. 

Given an arbitrary equation of state of the type p = p{n, s), the action functional is 



\J\ 



,s) + 



+ j^'id^^ + sd^e + A^M^i}- (A8) 

By varying the 4-metric we get the standard stress-energy-momentum tensor 

2 3'^ 

T"" = -i^ji = -epUi'U'' + p(V - dJi'U'') = -e(p + p)C/^C/'^ + p V^ (A9) 



where the pressure is given by 



dn 

The Euler-Lagrange equations for the fiuid motion are 

5S 



p = n^\s-p. (AlO) 



pU^ + d^if + sd^e + I3id^dt = 0, 
-d^J^ = 0, 
-9^(sJ^) = 0, 



5S_ 

5S 



Ji'd^a' = 0. (All) 



The second equation is the particle number conservation, the third one the entropy 
exchange constraint and the last one restricts the fiuid 4-velocity vector to be directed 
along the fiow lines a* = const.. The first equation gives the Clebsch or velocity-potential 
representation of the 4-velocity [/^ (the scalar fields in this representation are called Clebsch 
or velocity potentials). The fifth equations imply the constancy of the /9j's along the fiuid fiow 
lines, so that these Lagrange multipliers can be expressed as a function of the Lagrangian 
coordinates. The fourth equation, after a comparison with the first law of thermodynamics, 
leads to the identification T — U^^d^O = |n for the fiuid temperature. 

Moreover, one can show that the Euler-Lagrange equations imply the conservation of 
the stress-energy- momentum tensor T'^'^.^, — 0. This equations can be split in the projection 
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along the fluid flow lines and in the one orthogonal to them: 



i) The projection along the fluid flow lines plus the particle number conservation give 
U^T^^ -y = —^U^dfj_s = 0, which is verified due to the entropy exchange constraint. There- 
fore, the fluid flow is locally adiabatic, that is the entropy per particle along the fluid flow 
lines is conserved. 



ii) The projection orthogonal to the fluid flow lines gives the Euler equations, relating the 
fluid acceleration to the gradient of pressure 

(V - eU,U,)T''^,^ = -e(p + p)U,,,U'' - (5^ - eU,U^)d,p. (A12) 

By using p = n^\s — p, it is shown in Ref. IQ that these equations can be rewritten as 

2(/x[/[^);,]f/- = -e(<5; - [/^f/-)i|^|„9,s. (A13) 

The use of the entropy exchange constraint allows the rewrite the equations in the form 

2%;,]f/^ = T9^s, (A14) 

where V"^ = /if/^ is the Taub current (important for the description of circulation and 
vorticity), which can be identified with the 4-momentum per particle of a small amount of 
fluid to be injected in a larger sample of fluid without changing the total fluid volume or 
the entropy per particle. Now from the Euler-Lagrange we get 

2%;Hf/'^ = -2(5[^<^ + sdy,e + = {s^^^.^),,p'' = Td,s, (A15) 

and this result implies the validity of the Euler equations. 

In the non-relativistic limit {nU^).^^ = 0, T^'^.^, = become the particle number (or mass) 
conservation law, the entropy conservation law and the Euler-Newton equations. See Refs. 
23| , |30[| for the post-Newtonian approximation. 



We refer to Ref. |I| for the complete discussion. The previous action has the advantage 
on other actions that the canonical momenta conjugate to (f and 6 are the particle number 
density and entropy density seen by Eulerian observers at rest in space. The action evaluated 



on the solutions of the equations of motion is / d'^z^'^g{z)p{z). 

In Ref. [m there is a study of a special class of global Noether symmetries of this action 
associated with arbitrary functions F{a, j3i, s). It is shown that for each F there is a con- 
servation equation d^{FJ^) = and a Noether charge Q[F] = (Pa^n {el^U^)F{a, (3i, s) 
[S is a spacelike hypersurface with future pointing unit normal and with a 3-metric with 
determinant For F = 1 inside a volume V in S we get the conservation of particle 

number within a flow tube defined by the bundle of flow lines contained in the volume V. The 
factor el^U^ is the relativistic 'gamma factor' characterizing a boost from the Lagrangian 
observers with 4-velocity to the Eulerian observers with 4-velocity thus n{dpU^) is 
the particle number density as seen from the Eulerian observers. These symmetries describe 
the changes of Lagrangian coordinates a* and the fact that both the Lagrange multipliers 
and 9 are constant along each flow line (so that it is possible to transform any solution 



48 



to the fluid equations of motion into a solution with (p = 6 = on any given spacelike 
hypersurface) . 

However, the Hamiltonian formulation associated with this action is not trivial, because 
the many redundant variables present in it give rise to many first and second class constraints. 
In particular we get: 

1) second class constraints: 

A) TTjr ^0, r-n^^ 0; 

B) TTs ^ 0, sJ^ -TTe^O; 

C) TT^, ~ 0, PsJ^ - rCas ^ 0. 

2) first class constraints: iTjr ^ 0, so that the J^^s are gauge variables. 

Therefore the physical variables are the five pairs: ip, n^; 9, ttb; a\ iiar and one could study 
the associated canonical reduction. 

In Ref. [see its rich bibliography for the references] there is a systematic study of the 
action principles associated to the three types of equations of state present in the literature, 
first by using the Clebsch potentials and the associated Lagrange multipliers, then only in 
terms of the Lagrangian coordinates by inserting the solution of some of the Euler-Lagrange 
equations in the original action and eventually by adding surface terms. 

1) Equation of state p = p{n, s). One has the action 



S[n, U'',Lp,e,s,a'',pr;^gt,u] 



d^xJ^g p{n, s) — nU^{d^ip — 6d^s + Prdf^ 



(A16) 



If one knows s = s(a^) and = J^{a'^) = —\/^e^^P'^dyOi^dpa^daa^r]i2z{oi^)-, one can 
define S = S — J d'^xdf^[{(p + s6)J^], and one can show that it has the form 



S = S[a''] 



d xd^gp{ 



Ml 



(A17) 



2) Equation of state: p = p{fi, s) [V^ = piU^ Taub vector] 



S{p) = S(p) [V, p>, e, s, Pr] ^g^,u] = 

I — r dv / V^^ 

d^x^^g[p{p, s) - ^[\V\ - ^^{d^^ + sd,e + 



(A18) 



or by using one of its EL equations V^= — {d^(p + sd^O + Prd^a^) to eliminate one gets 
Schutz's action [p determined by /i^ = — V^V^] 



S{p)[^,0,s,a\(3r] g^,u\ = J d x^^gp{p,s) 
3) Equation of state a = a{n,T). The action is 



(A19) 



J\ 



g 



(A20) 



or 5'(a)[v5,6',s,a^/5^] = - / J|a(-J^, '^^d^i 



At the end of Ref. [0 there is the action for "isentropic" fluids and for their particular 
case of a "dust" (used in Ref. |^ as a reference fluid in canonical gravity). 
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The isentropic fluids have equation of state a(n, T) = — sT with s = const, (constant 
value of the entropy per particle). By introducing (f' = Lp + s9, the action can be written in 
the form 



d X 



or 



Q r -r. 4 

'^{isentropic)\f^ i Qjiv 



V 9 



\J\ 



9 



(A21) 



(A22) 



The dust has equation of state p{n) = fin, namely a{n, T) = fi — sT so that we get zero 
dp 



pressure p = n^ — p = 0. Again with (p = ip + s9 the action becomes 



d X 



(A23) 



or with = —-{dfj,(p' + (3^3 ^a^) [In Ref. M = p,n rest mass (energy) density and 
T = if'/p, Wr = -Pr, = a- = -d^T + Wrd^Z^] 



S[,,,,^ [T, Z^ M, Wr-, V] = -lj d'x^9iH {U, VU. - e) , (A24) 



or 



S{dust)[oi'^', 9fiiy] 



d'^xp\J\ 



(A25) 



In Ref. [0 there is a study of the action ( |A24| ) since the dust is used as a reference fluid 
in general relativity. At the Hamiltonian level one gets: 

i) 3 pairs of second class constraints [vr||,(r, 5) ~ 0, 7r^j,(r, a) — Wr{T, (r)'KT{r, a) ~ 0], which 
allow the elimination of Wrir, a) and vr||,(r, (t); 

ii) a pair of second class constraints [ 7rM(T, cr) ^ plus the secondary M(r, a) — 

= (r, a) ~ 0], which allow the elimination of M, ttm- 
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APPENDIX B: COVARIANT RELATIVISTIC THERMODYNAMICS OF 
EQUILIBRIUM AND NON-EQUILIBRIUM. 



In this Appendix we shall collect some results on relativistic fluids which are well known 
but scattered in the specialized literature. We shall use essentially Ref. 0, which has to be 
consulted for the relevant bibliography. See also Ref. . 



Firstly we remind some notions of covariant thermodynamics of equilibrium. 

Let us remember that given the stress-energy-momentum tensor of a continuous medium 
T'^'^, the densities of energy and momentum are T°° and c~^T^° respectively [so that dP^ = 
c'^rjT^'^dTju is the 4-momentum that crosses the 3-area element dHy in the sense of its 
normal [t] = —1 if the normal is spacelike, r] = +1 if it is timelike)]; instead, cT°^ is the 
energy flux in the positive r direction, while T^** is the r component of the stress in the plane 
perpendicular to the s direction (a pressure, if it is positive). A local observer with timelike 
4- velocity {v? = ec^) will measure energy density c~'^T^^^u^Uy and energy flux eT^^u^^riy 
along the direction of a unit vector in his rest frame. 

For a fluid at thermal equilibrium with T'^'' = pU^'U'' - e^{^g^''' - eU'^U'') [U'' is the 
hydrodynamical 4-velocity of the fluid] with particle number density n, specific volume 
V = - and entropy per particle s = (fcs is Boltzmann's constant) in its rest frame, the 
energy density is 

p(? = n{m(? + e), (Bl) 

where e is the mean internal (thermal plus chemical) energy per particle and m is particle's 
rest mass. 

From a non-relativistic point of view, by writing the equation of state in the form s = 
s{e, V) the temperature and the pressure emerge as partial derivatives from the first law of 
thermodynamics in the form (Gibbs equation) 

ds{e,V) = ^{de + pdV). (B2) 

If facias = e+pV — Ts is the non-relativistic chemical potential per particle, its relativistic 
version is 

p = mc^ + ficiass = p-Ts, (B3) 

[/i = is the specific enthalpy, also called chemical potential as in Appendix A] and we 

get 

/i n = pc^ + p — nTs = pc^ + p — ksTS^ 
ksTdS = d{pc^) — p dn = d{pc^) — {p — Ts)dn, or 
d{p(?) = p dn + Td{ns) = pdn + riTds. (B4) 

By introducing the "thermal potential" a = = and the inverse temperature 

2 

f3 = these two equations take the form 



51 



^ ns ^, p , 
S = — = (3[p + —)- an, 
Kb c 



dS = (3dp — adn. (B5) 

Let us remark that in Refs. [^,^,^] one uses different notations, some of which are 
given in the following equation (in Ref. p is denoted e and p is denoted r) 

p + 4 = n{mc^ + e) + ^ = p'h = p {c^ + e + = p {c^ + h'), (B6) 

where p = nm is the rest-mass density [r^: = y/^p is called the coordinate rest-mass 
density] and e = e/m is the specific internal energy [so that p h is the "effective inertial 
mass of the fluid; in the post-Newtonian approximation of Ref. p3| it is shown that a = 
c-^{T°" + T'') + 0(c-^) = c-2y4^(-T° + T/) + 0{c-^) has the interpretation of equality 
of the "passive" and the "active" gravitational mass] . For the specific enthalpy or chemical 
potential we get [p/m = h = + h is called enthalpy] 

p = -{p+P-) = ^{p +P^)=mh = m{c^ + h'). (B7) 
n p & 

See Ref. for a richer table of conversion of notations. 

Relativistically, we must consider, besides the stress-energy-momentum tensor T^^ and 
the associated 4-momentum = d^Y^^T^^ , a particle flux density (one nf^ for each 
constituent a of the system) and the entropy flux density s^. At thermal equilibrium all 
these a priori unrelated 4-vectors must all be parallel to the hydrodynamical 4-velocity 

= nf/^ , = sf/^, P'^ = PU^", (B8) 



2 

Analogously, we have V'^ = VU^ (V = 1/n is the specific volume), = (3U^ = j^U^ 

[a related 4- vector is the equilibrium parameter 4- vector = p (3^]. 

Since ell^T^'^ = pU", we get the final manifestly covariant form of the previous two 
equations (now the hydrodynamical 4-velocity is considered as an extra thermodynamical 
variable) 



dS" = -adn" - eMT"^- (B9) 
Global thermal equilibrium imposes 

d^a = d^(3, + d,(3^ = 0. (BIO) 

As a consequence we get 

di^p^") = n^'da + eT'^^dp^, (Bl 1) 

namely the basic variables n", T^" and S'^ can all be generated from partial derivatives of 
the "fugacity" 4- vector (or "thermodynamical potential") 
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^'=^-£^ ^(--) = |£' = r - c.n^ - ^.Ti::^at)^ (B12) 

once the equation of state is known. Here T^^^^) is the canonical or material (in general non 
symmetric) stress tensor, ensuring that reversible flows of field energy are not accompanied 
by an entropy flux. 

This final form remains valid (at least to first order in deviations) for "states that deviate 
from equilibrium", when the 4-vectors 5^, n^,... are no more parallel; the extra information 
in this equation is precisely the standard linear relation between entropy flux and heat flux. 
The second law of thermodynamics for relativistic systems is d^j^S^ > 0, which becomes a 
strict equality in equilibrium. 

The fugacity 4-vector is evaluated by using the covariant relativistic statistical theory 
for thermal equilibrium starting from a grand canonical ensemble with density matrix p 
by maximizing the entropy S = —Tr{pln p) subject to the constraints Trp = 1, Tr{ph) = n, 
Tr{pP'^) = P^: this gives (in the large volume limit) 

with 

lnZ=[ erdJ:^, n= I en^dT.^, P^' = f eTlZat)dK, (B13) 

JAS JAS JAT. ^ ' 

[it is assumed that the members of the ensemble are small (macroscopic) subregions of one 
extended body in thermal equilibrium, whose worldtubes intersect an arbitrary spacelike 
hypersurface in small 3-areas AE]. Therefore, one has to find the grand canonical partition 
function 

Z{V^,(3^,t^)=Y,e''^^'Q^{V„f3,), where Qn{V,,P,)=l rfa„(g,p)e-^''^", (B14) 

is the canonical partition function for fixed volume and dcTn is the invariant microcanonical 
density of states. For an ideal Boltzmann gas of N free particles of mass m [see Section III 
for its equation of state] it is 

I r N N 

danip, ^) = 1^ ^'(P -T.Pi)I[ ^V,P'0{pmPl - em')d'p,. (B15) 
■ •' 1=1 i=i 



Following Ref. [^ [using a certain type of gauge fixings to the first class constraints 
•pi — em^ ^ 0] in Ref. [^ Qn was evaluated in the rest-frame instant form on the Wigner 
hyperplane (this method can be extended to a gas of molecules, which are N-body bound 
states): 

Qn = J^[^K,{mP)] . (B16) 

The same results may be obtained by starting from the covariant relativistic kinetic 
theory of gas [see Ref. [^,|D|; in Ref. [^ there is a short review] whose particles interact 
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only by collisions by using Synge's invariant distribution function N{q,p) ||4T|] [the number of 
particle worldlines with momenta in the range {p^, duo) that cross a target 3-area dE^ in 
in the direction of its normal is given by dN = N{q,p)dujriv^dTi^ (= Nd^qd^p for the 3-space 
q° = const.); duj = d^p/v°^/^ is the invariant element of 3-area on the mass-shell]. One 
arrives at a transport equation for N,^ = V^{Nv'^) [v^ is the particle velocity obtained from 



the Hamilton equation implied by the one-particle Hamiltonian H = y €'^g^^'^{q)pf^Pi, = m 
(it is the energy after the gauge fixing g° r to the first class contraint '^g^"{q)PnPv — ~ 
0); is the covariant gradient holding the 4- vector p^ (not its components) fixed] with 
a "collision term" C[N] describing the collisions; for a dilute simple gas dominated by 
binary collisions one arrives at the Boltzmann equation [for C[N] = one solution is the 
relativistic version of the Maxwell-Boltzmann distribution function, i.e. the classical Jiittner- 
Synge one N = const.e'^^'^^ / 'i'Km'^K2{m(3) for the Boltzmann gas |^1|]]. The H-theorem 
[V^S*^ > 0, where S^{q) = — J[Nln{Nh^) — N]v^duj is the entropy flux] and the results 
at thermal equilibrium emerge [from the balance law V^(/ N fv^duj) = J fC[N]duj (/ is 
an arbitrary tensorial function) one can deduce the conservation laws V^n^ = ^T'^^ = 0, 
where n'^ = J Nv^duj, T^^ = J Npi,v^duj; the vanishing of entropy production at local 
thermal equilibrium gives Neq{q,p) = /;,~3ga(g)+/3i,(g)P"' jj-^ ^^le case of Boltzmann statistic 
and one gets (?7^ = /3^'//3) n^^ = J NegVi'diu = nU^", T^^ = pU''U>' - epi^g'''' - eWW), 
^eq = vP^ ~ ct'^eq ~ PuT^q ouc obtaius the equations for the Boltzmann ideal gas given in 
Section III]. 

One can study the small deviations from thermal equilibrium [A^ = Neq{l + f), where Neq 
is an arbitrary local equilibrium distribution] with the linearized Boltzmann equation and 
then by using either the Chapman- Enskog ansatz of quasi-stationarity of small deviations 
(this ignores the gradients of / and gives the standard Landau-Lifshitz and Eckart phe- 
nomenological laws; one gets Fourier equation for heat conduction and the Navier-Stokes 
equation for the bulk and shear stresses; however one has parabolic and not hyperbolic 
equations implying non-causal propagation) or with the Grad method in the 14-moment 
approximation. This method retains the gradients of / [there are 5 extra thermodynami- 
cal variables, which can be explicitly determined from 14 moments among the infinite set 
of moments / Np'^p'^p''...d'^p of kinetic theory; no extra auxiliary state variables are intro- 
duced to specify a non-equilibrium state besides T^'^, n^, S^] and gives phenomenological 
laws which are the kinetic equivalent of Miiller extended thermodynamics and its various 
developments; now the equations are hyperbolic, there is no causality problem but there 
are problems with shock waves. See Ref. [^ for the bibliography and for a review of the 
non-equilibrium phenomenological laws (see also Ref. |^2[) of Eckart, Landau-Lifshitz, of the 
various formulations of extended thermodynamics, of non-local thermodynamics. 

While in Ref. |^3[ it is said that the difference between causal hyperbolic theories and 
acausual parabolic one is unobservable, in Ref. [^[see also Ref. [^] there is a discussion 
of the cases in which hyperbolic theories are relevant. See also the numerical codes of Refs. 



mm 

In phenomenological theories the starting point are the equations d^T'^'^ = dfj_n^ = 0, 
d^S'^ > 0. There is the problem of how to define a 4-velocity and a rest-frame for a given 
non-equilibrium state. Another problem is how to specify a non-equilibrium state completely 
at the macroscopic level: a priori one could need an infinite number of auxiliary quantities 
(vanishing at equilibrium) and an equation of state depending on them. The basic postulate 
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of extended thermodynamics is the absence of such variables. 

Regarding the rest frame problem there are two main solutions in the literature connected 
with the relativistic description of "heat flow": 

i) Eckart theory. One considers a local observer in a simple fluid who is at rest with 
respect to the average motion of the particles: its 4-velocity ?7(ecfc) is parallel by definition 
to the particle fiux n^, namely 

= n(ecfc)t/fec.). (B17) 

This local observer sees "heat fiow" as a fiux of energy in his rest frame: 
SO that we get 

-t — P{eck)'-^leck)'-^{eck) + 'i{eck)'-^ (eck) ^" {eck)'i{eck) + ^{eck)-' 
PteO^) = Pie'ck) = <P + neck)){V' ' + 

P(Zk)Uieck)i^ = (l{eck),jU'(^eck) = 0' '^{eck)nv{^ 9^"" - eU'^^^f.-^Ul^^j,)) = 0, (B18) 

where p is the thermodynamic pressure, 7r(^eck) the bulk viscosity and tt^^J^^^ the shear stress. 
This description has the particle conservation law d/^n'^ = 0. 

ii) Landau- Lifshitz theory. One considers a different observer (drifting slowly in the 
direction of heat fiow with a 3- velocity vd — q/nmc?) whose 4-velocity Uq^i-^ is by definition 
such to give a vanishing "heat fiow", i.e. there is no net energy flux in his rest frame: 
Ul^ii^Tl^nu = for all vectors orthogonal to U^^y This implies that Uj^n-^ is the timelike 
eigenvector of T^*^, T^'^Uqi^i, = epQi^Uj^ny which is unique if T'^'^ satisfies a positive energy 
condition. Now we get 

n" = n^ii)U^ii) + j^ii), J(ii)f.U^ii) = (J = -nvD = -q/mc^). (B19) 

This observer in his rest frame does not see a heat fiow but a particle drift. This 
description has the simplest form of the energy-momentum tensor. 

One has n^^ck) = n^ii)ch(p, pi^^ck) = p{ii)ch^ ^ + P{u)sh^ V = t^^" 3ti3v/n%ckY with chip ^ 

^(ii)^ieck)^i [the difference is a Lorcntz factor \Jl — vj^/c^, so that there are insignificant 
differences for many practical purposes if deviations from equilibrium are small]. The angle 
if K. jjn ~ vd/c ~ qjnm(? is a dimensionless measure of the deviation from equilibrium 
{ni^ii) - n^eck) and p^u) - p(^^ck) are of order cp'^]. 
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One can decompose T^^, n^, in terms of any 4-velocity that falls within a cone of 
angle ^ (f containing Uj^^^i^^ and Ufl^y, each choice gives a particle density n{U) = eu^n'^ 
and energy density p{U) = U^U^T'^'^ which are independent of if one neglects terms of 
order cp'^. Therefore, one has: 

i) if Seq{p{U),n{U)) is the equilibrium entropy density, then S{U) — eU^S^ — Seq + 0{(p'^); 

ii) if p{U) = — l^/^ls/n is the (reversible) thcrmodynamical pressure defined as work done 
in an isentropic expansion (off equilibrium this definition allows to sepate it from the bulk 
stress 7i{U) in the stress-energy-momentum tensor) and Peq is the pressure at equilibrium, 
then p{U) = Peq{p{U),n{U)) + 0{ip''). 

By postulating that the covariant Gibbs relation remains valid for arbitrary infinitesimal 
displacements (Sn'^.dT^'^ , ..) from an equilibrium state, one gets a covariant off-equilibrium 
thermodynamics based on the equation 

= p{a, 0)0" - an" - fi^T"" - Q"{5n\ ST"", ..), 
V^^S" = -5n"d^a - 5T""V.P^ - W^Q" > 0, (B20) 

with Q" of second order in the displacements and a, (5^ arbitrary. At equilibrium one recovers 
^ _ _ I3^T^^, V^S^g = [with U" = P"/P and (for viscous heat-conducting 

fluids, but not for superfiuids) d^a = V^Pu + = 0] . 

If we choose P" = U" /ksT parallel to n" of the given off-equilibrium state, we are in the 

"Eckart frame" , U" = U"^^^^ , and we get 

5" = €U(^eck)ixS" = Seq + dJ{eck)nQ^ , 
ieck) = -{"r - eU^eok)UU))nU^eck),. (B21) 

SO that to linear order we get the standard relation between entropy flux df^^ck) and heat flux 

?(ecfe) 

^(ecfe) = + {possible 2nd order term). (B22) 

If we choose U" = U^u), the timclike eigenvector of T"" , so that UQi^^Tl^C^gp — 
€-Uqi)U{ii)p) — 0, we are in the "Landau- Lifshitz frame" and we get 

j'(u)-{'r-eU^u)Ula))n.. (B23) 
so that at linear order we get the standard relation between entropy flux and diffusive 

flux Jill) 
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= —-J — + {possible 2nd order term). (B24) 

In the Landau-Lifschitz frame heat flow and diffusion are englobed in the diffusive flux 
Jqi) relative to the mean mass-energy flow. 



The entropy inequality becomes (each term is of second order in the deviations from lo- 
cal equilibrium) 

< V^^'^ = -Sn'^di.a - ST'^^V^Pi, - V^Q^ (B25) 

with the fitting conditions Sn^Ufj, — ST^^Uf^U,, — 0, which contain all information about the 
viscous stresses, heat fiow and diffusion in the off-equilibrium state (they are dependent on 
the arbitrary choice of the 4- velocity U^). 

Once a detailed form of is specified, linear relations between irreversible fiuxes ST'^", 
Sn^^ and gradients V(^/5,^), d^a follow. 

A) Q^^ — (hke in the non-relativistic case) . 

The spatial entropy fiux a is only a strictly linear function of heat flux q and diffusion 
flux j. In this case the off-equilibrium entropy density S = eU^S^ is given by the equilib- 
rium equation of state 5" = Seq{p,n). We have < V^S''* = —Sn'^df^a — 5T'^''Vj,(5^ with 
fltting conditions Sn^U^ — ST^^'^UnU„ — and with still arbitrary at flrst order. 

Al) Landau-Lifschitz frame and theory. If^ = U^^-^ is the timclikc eigenvector of T^'^. 
This and the fltting conditions imply 5T^^U(ii)u = 0. The shear and bulk stresses n'^^^ tt^u) 
are identifled by the decomposition 

The inequality V^5''' > becomes 

- \{'9,. - eU^u),U^^u)u){'g-^ - eU^u)U^ii^)]X^p, (B27) 
[the < .. > operation extracts the purely spatial, trace- free part of any tensor]. 

If the equilibrium state is isotropic (Curie's principle) and if we assume that {3q,i)iT^qj)iT^{ii)) 
are "linear and purely local" functions of the gradiants, V^S*^ > implies 

3'ii)--<9"'-eU^u)U'(u))d,a, > 0, (B28) 

[it is a mixture of Fourier's law of heat conduction and of Pick's law of diffusion, stemming 
from the relativistic mass-energy equivalence]. 
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and the standard Navier-Stokes equations {(s, Cv are shear and bulk viscosities) 

vr(H)M- = -2Cs < V./3^ >, 71(11) = ^CyV^t/J,). (B29) 

A2) Eckart frame and theory. U^^^,^-^ parellel to n'*. Now we have the fitting condition 
Sn'^ = 0. The heat fiux appears in the decomposition of ST^'^ [0'{eck)n — U^eck)^ i'^{eck)n is 
the 4-acceleration] 

ST'-" = qU)UU) + Kck)^(eck) + <ck) + ^^eck){V'' ' ^eckfU)) ■ (B30) 

The inequality V^^^ > becomes 

<ecfe)(^M« - /5a(ecik)M) " /^(Trfecfe) V.C/(ecfe)M + ^{eck)^ ^^%ck)) ^ 0" (^31) 

With the simplest assumption of linearity and locality, we obtain Fourier's law of heat 
conduction [it is not strictly equivalent to the Landau-Lifshitz one, because they differ by 
spatial gradients of the viscous stresses and the time-derivative of the heat fiux] 

qU) = -^'a'" - eU^eck)U;:eck)){d.T + mC/(ee.).), (B32) 

[the term depending on the acceleration is sometimes referred to as an effect of the "inertia 
of heat"], and the same form of the Navier-Stokes equations for n'^^^f^^, '^{eck) (they are not 
strictly equivalent to the Landau-Lifshitz ones, because they differ by gradients of the drift 
vd — q/nmc^). 

For a simple fluid Fourier's law and Navier-Stokes equations (9 equations) and the con- 
servation laws V ^T^^ = Vj^n'* = (5 equations) determine the 14 variables T'^'' , from 
suitable initial data. However, these equations are of mixed parabolic-hyperbolic- elliptic 
type and, as said, one gets acausality and instability. 

Kinetic theory gives 

Q'--\j NejVdu; ^ 0, (B33) 

for a gas up to second order in the deviation (A^ — N^q) = N^qf [Q'^ ~ requires small gra- 
dients and quasi-stationary processes] . Two alternative classes of phenomenological theories 
are 

B) Linear non-local thermodynamics (NLT). 

This theory gives a rheomorphic rather than causal description of the phenomenological 
laws: transport coefficients at an event x arc taken to depend, not on the entire causal 
past of X, but only on the past history of a "comoving local fluid element". It is a lin- 
ear theory restricted to small deviations from equilibrium, which can be derived from the 
linearized Boltzmann equation by projector-operator techniques (and probably inherits its 
causality properties). Instead of writing ((5n^(a;), (5T^j,(x)) = (^{U, T){—dfj_a{x), — V(^/9i,)(a;)), 
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this local phenomenological law is generalized to ((5n^(x, 5T^jy(x, = J^^dx°' a{x° 



C) Local non-linear extended thermodynamics (ET). 

It is more relevant for relativistic astrophysics, where correlation and memory effects are 
not of primary interest and, instead, one needs a tractable and consistent transport the- 
ory coextensive at the macroscopic level with Boltzmann's equation. It is assumed that the 
second order term Q'^{6n^, 5T"p, ..) does not depend on auxiliary variables vanishing at equi- 
librium: this ansatz is the phenomenological equivalent of Grad's 14-moment approximation 
in kinetic theory. These theories are called "second-order theories" and many of them are 
analyzed in Ref. |^6|]; when the dissipative fluxes are subject to a conservation equation. 



these theories are called of causal "divergence type" like the ones of Refs. ||47|-^ . Another 
type of theory (extended irreversible thermodynamics; in general these theories are not of 
divergence type) was developed in Refs. ||50H53|: in it there are transport equations for the 
dissipative fluxes rather than conservation laws. 

For small deviations one retains only the quadratic terms in the Taylor expansion of 
(leading to "linear" phenomenological laws): this implies 5 new undetermined coefficients 

= \u^[l3oT^^ + I3iqy, + /327r^'^7r^.] - ^ovrg^ - aivr^'^^g., (B34) 

with /3i > from eU^Q'^ > (the /3i's are 'relaxation times'). A first-order change of 
rest frame produces a second-order change in Q'^ [going from the Landau-Lifshitz frame 
to the Eckart one, one gets a(eck)i - Oi{ii)i = P{ii)i - P{eck)i = [(p + p)T]-^ /3(zz)o = P(eck)o, 
P{ii)2 = /3(ecfe)2, and the phenomenological laws are now invariant to first order]. 

In the "Eckart frame" the phenomenological laws take the form 

(l(eck) = -I^T{^g^'' - et/(lfc)f/(ecfc))[r~^<9^^ + a(e,fc)^ + P(eck)ldrq{eck)u " 
— 0![eck)odu'TC{eck) — Q^(ecfc)l VpTT^'g^^^j^] , 

'^{eck)^iu = —2Cs[< ^uU(eck)^i > + P {eck)2dTTC {eck) ^lu — Oi(eck)l < Vyg(ecfc)/^ >], 

7r(ecfc) = -^Cvl"^ fiUl^eck) + P(eck)odr7l(eck) " «(ecfc)o V^gf^^^)] , (B35) 

which reduce to the equation of the standard Eckart theory if the 5 relaxation (/3j) and 



coupling (ctj) coefficients are put equal to zero. See for instance Ref. ||5^ for a complete 
treatment and also Ref. [^]. For appropriate values of these coefficients these equations are 
hyperbolic and, therefore, causal and stable. The transport equations can be understood 



as evolution equations for the dissipative variables as they describe how these fluxes 



evolve from an initial arbitrary state to a flnal steady one [the time parameter r is usually 
interpreted as the relaxation time of the dissipative processes] . In the case of a gas the new 
coefficients can be found explicitly [see also Ref. ||56| for a recent approach to relativistic 
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interacting gases starting from the Boltzmann equation], and they are purely thermody- 
namical functions. Wave front speeds are finite and comparable with the speed of sound. 
A problem with these theories is that they do not admit a regular shock structure (like the 
Navier-Stokes equations) once the speed of the shock front exceeds the highest characteristic 
velocity (a "subshock" will form within a shock layer for speeds exceeding the wave-front 
velocities of thermo- viscous effects). The situation slowly ameliorates if more moments are 
taken into account |B3 . 



In the approach reviewed in Ref. the extra indeterminacy associated to the new 5 
coefficients is eliminated (at the price of high non-linearity) by annexing to the usual con- 
servation and entropy laws a new phenomenological assumption (in this way one obtains a 
causal divergence type theory): 

^^A^^p = jA*^^ (B36) 

in which A^'^'^ and I'^" are symmetric tensors with the following traces 

A^^^ = -n", J'^^ = 0. (B37) 

These conditions are modelled on kinetic theory, in which A'''^" represents the third mo- 
ment of the distribution function in momentum space, and I'^'^ the second moment of the 
collision term in Boltzmann's equation. The previous equations are central in the determina- 
tion of the distribution function in Grad's 14-moment approximation. The phenomenologi- 
cal theory is completed by the postulate that the state variables S'^, A^^'^ , I'^'^ are invariant 
functions of T^*^, only. The theory is an almost exact phenomenological counterpart of 
the Grad approximation. See Ref. ||3^ for the beginning (only non viscous heta conduct- 
ing materials are treated) of a derivation of extended thermodynamics from a variational 
principle. 

Everything may be rephrased in terms of the Lagrangian coordinates of the fluid used in 
this paper. What is lacking in the non-dissipative case of heat conduction is the functional 
form of the off-equilibrium equation of state reducing to p = p{n, s) at thermal equilibrium. 
In the dissipative case the system is open and T'^'^, P^, are not conserved. 

See Ref. for attempts to deflne a classical theory of dissipation in the Hamiltonian 



framework and Ref. about Hamiltonian molecular dynamics for the addition of an extra 
degree of freedom to an N-body system to transform it into an open system (with the choice 
of a suitable potential for the extra variable the equilibrium distribution function of the 
N-body subsystem is exactly the canonical ensemble). 

However, the most constructive procedure is to get (starting from an action principle) 
the Hamiltonian form of the energy-momentum of a closed system, like it has been done in 



Ref. |I3 for a system of N charged scalar particles, in which the mutual act ion- at- a- distance 
interaction is the complete Darwin potential extracted from the Lienard-Wiechert solution in 
the radiation gauge (the interactions are momentum- and, therefore, velocity-dependent). In 
this case one can deflne an open (in general dissipative) subsystem by considering a cluster of 
n < N particles and assigning to it a non-conserved energy-momentum tensor built with all 
the terms of the original energy-momentum tensor which depend on the canonical variables 
of the n particles (the other N — n particles are considered as external flelds). 
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APPENDIX C: NOTATIONS ON SPACELIKE HYPERSURFACES. 



Let us first review some preliminary results from Refs. [|T^ needed in tlie description of 
pliysical systems on spacelike liyper surf aces. 

Let {S^} be a one-parameter family of spacelike hypersurfaces foliating Minkowski space- 
time M*^ with 4-metric rj^^ = e(H ), e = ± [e = +1 is the particle physics convention; 

e = — 1 the general relativity one] and giving a 3-M decomposition of it. At fixed r, let 
z'^{t, a) be the coordinates of the points on E,- in M^, {a} a system of coordinates on S^. 
If = (o""^ = T;cr = {a^}) [the notation A = (r, f) with f = 1, 2, 3 will be used; note that 
A = T and A = f = 1,2,3 are Lorentz-scalar indices] and = d/da"^, one can define the 
vierbeins 

4(r, a) = d^z^ir, a), d^z\ - d^z% = 0, (CI) 

so that the metric on St- is 

dABi.T^ ^) = Zj^i^, ^)Vf^uz'"^{T, a), egrrir, a) > 0, 
g{T,a) = -det\\g^BiT,a) \ \ = {det\\ z')^{t, a) \ 

'j{T,a) = -det\\gfsiT,a) \\ = det\\^gfsir, a)\\, (C2) 

where gfs = — e^fi'rs with ^gfs having positive signature (+ + +). 

If Y'^ijj ^) = —^^g^'^ is the inverse of the 3-metric gfs{T, a) [7™(t, (T)gus{j, B) = 6^], the 
inverse g^^{T, a) of gAsiT, ^) ^)9cb(j, ^) = ^i] is given by 



g^'{T,a) = -[2^,,7^^](r,a) = e[^g,^' g'^%T,a), 
9 9 

/«-(r, a) = Y\r, a) + [lg^^g^,^^^'^f^]{r, a) = 

-e Y'(r,a) + [-gr^gr,'9uf'g'']iT, a), (C3) 
9 



so that 1 = g {t, a)g(j^{T, a) is equivalent to 



We have 



and 



% = 9Tr{T, a) - (r, a) g^f (r, a) g^s (r, a) . ( C4) 

7(r, a) 



^'Ar,a) = r^l^ + g^,Y'z^){r,a), (C5) 
V 7 



= (Pr + z,V%^)(r,a), (C6) 
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where 

/^(r,a^) = (-^eW^f4)(r,a), 

/2(r,a) = l, /^(r,a)znT,f?) = 0, (C7) 

is the unit (future pointing) normal to at z'^{t, a). 

For the volume element in Minkowski spacetime we have 

d'^z = z^ij, a)dT<f'L^ = dT[z'^{T, 3)1^,{t, (y)\\J ^{T,a)d^a = 

g{T, a)dTd^a. (C8) 



Let us remark that according to the geometrical approach of Ref. [0,one can use Eq.(|C^ 
in the form 

z^f (r, a) = iV(r, a)/^(r, a) + W\t, a)z,^(r, a), 



where N = ^Jg/j = grr - Y^drfgrs = V9tt + e^g^^grrgrs and = grsY^ = -egrs ^g^^ 
are the standard lapse and shift functions N[z](fiat), N[z\{fiat) Introduction, so that 



g^r = eN^ + grsN'N' = e[N' - ^gfsN'N% 
Qrf = 9fsN' = -e^gfsN', 
g-- = eN-\ 
g^f = -eN'-/N^, 



d _ J _d_ _|_ ^ ^ sf d _ J _d_ 3„sr_9_ 



d^z = N^drd^a. 

The rest frame form of a timelike fourvector isP ^ = rjy/ep'^ll] 0) = ri'^°rj^/ep^, P"^ = p'^, 

o 

where rj = signp°. The standard Wigner boost transforming P into p'^ is 



l''Ap.p) = ^'uHp)) 



ep p. p +ep2 



1 + u°{p) 



,2 



= e^^(n(p)) = u^{p) = p^ /r]\Jep 
y = r e^Hp)) = {-Ur{p);6l-^p^^). (C9) 

o o 

The inverse of L^y{p,p) is L^y{p,p), the standard boost to the rest frame, defined by 
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L^(p,p) = W{p,p) = L^(p,p)|p-_^_^-. (CIO) 

Therefore, we can define the following vierbeins [the e^('u(p))'s are also called polarization 
vectors; the indices r, s will be used for A=l,2,3 and o for A = o] 

4i<P)) = L\{p,p) = 77^^77,.6^(^.(p)), 
^U'^ip)) = ViJ.vel{u{p)) = u^{p), 

elHp)) = -5^VAu{p)) = {6-Us{p); 5} - ^''^jHy^^^), 
eiiu{p))^uM, (Cll) 



which satisfy 



4{u{p)yBiu{p)) ^ vi 



r=l 



Vab = ^A{u{p))VnueB{u{p)), 

Pa^e'Xiuip)) = Pa^e^Hp)) = 0. (C12) 

dpa OPa 



The Wigner rotation corresponding to the Lorentz transformation A is 

i(A,p) i ■ 



B'„(A,p) = lL(P,p)A-'L{Ap,P)]\ = J ^.^ " 



Pp(A-i)''o + r/v^ 



p° + r^y/ep^^^ Ppik-^Yo + vV^ 

The polarization vectors transform under the Poincare transformations (a, A) in the 
following way 

e(f(i.(Ap)) = (i?-^)/A^er(ii(p)). (C14) 
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APPENDIX D: MORE ON DIXON'S MULTIPOLES. 



Let us add other forms of the Dixon multipoles. 

In the case of the fluid configurations treated in Section II and IV, the Hamilton equations 
generated by the Dirac Hamiltonian (|2.49| ) in the gauge g^ys ~ [A(r) = 0] imply [in Ref. 
20| this is a consequence of d^T^^'' = 0] 



Let us define for n > 1 



0, for n = 0, 



(Dl) 



= e(^n^(Ps))....6^f:(n(p.))e^(n(p,))/---^^(T.), 



n + 1 



^^•■■^"'^(T,) = cfi-''"^''(T,) =p^^-'^"'^(T,) -pV^^-^"''^(T,) = 
-e(f(«(p,))...6(f:(n(p,))e^(n(p,))]/---^^(T,), 



n 



(T.) = 0, 



n + 1 



(D2) 



and then for n> 2 



+ 



n 



n 



{u{Ps)) 



^{^ll...^ln^l)v 



{Ts) = 0, 
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T) — 1 

n 



+ e'sAu{ps))e:MPs)WT (Ts)- 



_d_ 



n 

+ 4^i-""^i^^)(T,)]. (D3) 
Then Eqs. (pl|) may be rewritten in the form 
1) n = l 

(T,) = p^^{T,)u^\p,) + = P^u^{ps)u^{ps) + P'u^^{ps)e';\u{p,)) + 

+ ei^{u{p.W\psW^{Ts) + 
+ ^Hn{Ps)VXn{Ps))lT^{Ts), 



2) n = 2 [identity fr,^'' = 1^^^^" + 4'''^^ + t^^''^<'^ 



T 

3) n > 3 



■^'■■■^-^^{T,)u''\p,) + -3^[^^6^^-^"^'(T,) + -4^-^"(^")(T,)], (D4) 



n dig n + 1 n 

This allows pO] to rewrite < T^", / > in the following form 



(2vr)^ 

n 

+ T.'-ij-K.---KJ^'---'"'''iTs)], (D5) 



n=2 
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with 



q-tl\...flntll' 

_L,fj-i 



(/il.../i„)(/il') 



n 



n 



n + 1 



n 



+ e. 



2 



n + 2 



^0^1. ..MnM)^^^^^ = 0, 



Ml 



n + 3 

n + 1 
1 



n 

+ <(«(p.))<.(«b.))[/T^-''"""(7;)- 



n 



. T(r\...rnS\S2) 



{Ts)]. 



Finally, a set of multipoles equivalent to the Xj}'"'^''^'' is 

n > 



n+1 



u 



l^l^l.-./ln [PO"] 



+ u^'{Ps)p'?'''-^"^^"\Ts) 



1 



(D6) 



^ ^ ^ [u^^{ps)e^Xu{Ps))e['{u{ps))e7{u{Ps)) + 
+ K[^(p,)6,-l(K(p,))e[''(i.(p,))e:i(K(p,)); 

[(n + 4)(3n + 5) linearly independent components], 



n>l 
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n>2 



4(n - I] 



J. 



n 



T 



"Ml " 



ri..rnAB 



n + l 

+ u^^{ps)e'^Xu{Ps))e^t{u{ps)VKu{ps)) 



Trri..rnsAT 



(Ts). 
(D7) 



The J^^ 



.pnpVpCT 



are the Dixon "2"+^-pole inertial moment tensors" of the extended sys- 



tem: they [or equivalently the Xj}"'^""^^^s\ determine its energy-momentum tensor together 
with the monopole pj. and the spin dipole Sq^ . The equations d^T^^^ = are satisfied due to 
the equations of motion ( p4| ) for and Sj^ [the so called Papapetrou-Dixon-Souriau equa- 
tions given in Eqs. (|4.18| )] without the need of the equations of motion for the ji^^---f^'^t^'^p'^ _ 
When all the multipoles jil:,i---i^'^i^'^p'^ ^re zero [or negligible] one speaks of a pole-dipole field 
configuration of the perfect fiuid. 
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